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SUMMARY 



In this thesis, we investigate the quantum mechanical interaction of hght with 
matter in the form of a gas of ultracold atoms: the atomic ensemble. We present 
a theoretical analysis of two problems, which involve the interaction of quantized 
electromagnetic fields (called signal and idler) with the atomic ensemble (i) cascade 
two-photon emission in an atomic ladder configuration, and (ii) photon frequency 
conversion in an atomic diamond configuration. The motivation of these studies 
comes from potential applications in long-distance quantum communication where it 
is desirable to generate quantum correlations between telecommunication wavelength 
light fields and ground level atomic coherences. In the two systems of interest, the 
light field produced in the upper arm of an atomic Rb level scheme is chosen to lie in 
the telecom window. The other field, resonant on a ground level transition, is in the 
near-infrared region of the spectrum. Telecom light is useful as it minimizes losses in 
the optical fiber transmission links of any two long-distance quantum communication 
device. 

We develop a theory of correlated signal-idler pair correlation. The analysis is 
complicated by the possible generation of multiple excitations in the atomic ensemble. 
An analytical treatment is given in the limit of a single excitation assuming adiabatic 
laser excitations. The analysis predicts superradiant timescales in the idler emission 
in agreement with experimental observation. To relax the restriction of a single 
excitation, we develop a different theory of cascade emission, which is solved by 
numerical simulation of classical stochastic differential equation using the theory of 
open quantum systems. The simulations are in good qualitative agreement with 
the analytical theory of superradiant timescales. We further analyze the feasibility 

xxi 



of this two-photon source to reahze the DLCZ protocol of the quantum repeater 
communication system. 

We provide a quantum theory of near-infrared to telecom wavelength conversion in 
the diamond configuration. The system provides a crucial part of a quantum-repeater 
memory element, which enables a "stored" near- infrared photon to be converted to a 
telecom wavelength for transmission without the destruction of light-atom quantum 
correlation. We calculate the theoretical conversion efficiency, analyzing the role of 
optical depth of the ensemble, pulse length, and quantum fluctuations on the process. 
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CHAPTER I 
INTRODUCTION 

A quantum communication network based on the distribution and sharing of entan- 
gled states is potentially secure to eavesdropping and is therefore of great practical 
interest [1, 2, 3]. A protocol for the realization of such a long distance system, 
known as the quantum repeater, was proposed by Briegel et al. [4, 5]. A quantum 
repeater based on the use of atomic ensembles as memory elements, distributed over 
the network, was subsequently suggested by Duan, Lukin, Cirac and Zoller [6]. The 
storage of information in the atomic ensembles involves the Raman scattering of an 
incident light beam from ground state atoms with the emission of a signal photon. 
The photon is correlated with the creation of a phased, ground-state, coherent ex- 
citation of the atomic ensemble. The information may be retrieved by a reverse 
Raman scattering process, sending the excitation back to the initial atomic ground 
state and generating an idler photon directionally correlated with the signal photon 
[7, 8, 9, 10, 11, 12, 13, 14, 15]. In the alkali gases, the signal and the idler field wave- 
lengths are in the near-infrared spectral region. This presents a wavelength mismatch 
with telecommunication wavelength optical fiber, which has a transmission window 
at longer wavelengths (1.1-1.6 um). It is this mismatch that motivates the search for 
alternative processes that can generate telecom wavelength photons correlated with 
atomic spin waves [16]. 

This motivates the research presented in this thesis where we study multi-level 
atomic schemes in which the transition between the excited states is resonant with a 
telecom wavelength light field [16]. The basic problem is to harness the absorption 
and the emission of telecom photons while preserving quantum correlations between 
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the atoms, which store information and the photons that carry along the optical fiber 
channel of the network. In this thesis, we theoretically study atomic cascade and 
diamond configurations in this context. 

1.1 DLCZ Protocol for the Quantum Repeater 

A long-distance quantum repeater must overcome the exponential losses in the optical 
fiber. To overcome this problem, the use of quantum memory was proposed [6]. For 
a practical system, it is essential to maximize quantum memory time, to preserve 
coherence during protocol operations, and connect the memory elements by light 
signals in the low-loss window of the optical fiber medium. The telecom wavelength 
range (1.1-1.6 fim) has a loss rate as low as 0.2 dB/km. 

It is not common to have a telecom ground state transition in atomic gases except 
for rare earth elements [17, 18] or in an erbium-doped crystal [19]. However, a 
telecom wavelength (signal) can be generated from transitions between excited levels 
in the alkali metals [16]. 

1.1.1 Correlated cascade emission in quantum telecommunication 

The ladder configuration of atomic levels provides a source for telecom photons (sig- 
nal) from the upper atomic transition. For rubidium and cesium atoms, the signal 
field has the range around 1.3-1.5 /im that can be coupled to an optical fiber and 
transmitted to a remote location. Cascade emission may result in pairs of photons, 
the signal entangled with the subsequently emitted infrared photon (idler) from the 
lower atomic transition. Entangled signal and idler photons were generated from a 
phase-matched four- wave mixing configuration in a cold, optically thick ^^Rb ensem- 
ble [16]. This correlated two-photon source is potentially useful as the signal field 
has telecom wavelength. 

The temporal emission characteristics of the idler field, generated on the lower 
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arm of the cascade transition, were observed in measurements of the joint signal- 
idler correlation function. The idler decay time was shorter than the natural atomic 
decay time and dependent on optical thickness in a way reminiscent of superradiance 
[23, 20, 21, 22, 24]. 

We will develop an analytical theory of the cascade emission in an atomic ensemble 
in Chapter 3. The influence of electromagnetic dipole-dipole interactions between 
atoms is important to account for the idler field's temporal profile. By developing the 
theory on the assumption of weak adiabatic laser excitation, we are able to calculate 
the spectral characteristics of the signal and idler fields, and make a connection with 
the traditional theory of superradiance. 

In Chapter 4, , we develop a more elaborate theory of the cascade emission under 
similar physical conditions to Chapter 3, but without the assumption of single atomic 
excitations. The theory is based on numerical solutions of stochastic differential 
equations derived using open-systems methods of quantum optics. We limit our 
analysis to the confirmation of the superradiant emission of the idler field predicted 
in the simple theory and observed experimentally. 

In Chapter 5, we use this theory to discuss a potential application of the cas- 
cade emission process in the DLCZ protocol, and discuss the role of time-frequency 
entanglement. 

1.2 Quantum Memory with Light Frequency Con- 
version 

It is not sufficient to generate telecom wavelength light for quantum communication. 
The light field must be quantum correlated with atomic excitations stored in memory 
[16]. 

Recently there has been a breakthrough in this direction using a pair of cold, non- 
degenerate rubidium gas samples [25]. A correlated pair of atomic spin wave and 
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infrared fields are generated by conventional Raman scattering in one ensemble. The 
light field is directed onto a second ensemble where it is frequency converted to the 
telecom range by four-wave mixing using a diamond configuration of atomic levels. 
The experiments were designed to measure quantum correlations between the stored 
atomic excitation and the telecom field. 

The conversion scheme exploits an efficient low-noise parametric conversion pro- 
cess that is facilitated by operating in the regime of high transparency [25]. This 
provides a basic quantum memory element for a scalable, long distance quantum net- 
work. In Chapter 6, we investigate conditions required to maximize the conversion 
efficiency as a function of optical thickness of the atomic ensemble. The infiuence of 
the probe pulse duration on the conversion efficiency is studied by numerical solution 
of the Maxwell-Bloch equations. 

1 . 3 Outline 

The remainder of this thesis is organized as follows. 

In Chapter 2, we review some theoretical methods to provide background for the 
theories developed in Chapter 4 and 6. In particular we discuss the derivation of 
quantum Heisenberg-Langevin equations for the interaction of a group of atoms with 
a quantized propagating electromagnetic field. We illustrate the connection of these 
operator equations with related classical (c-number) stochastic Langevin equations. 
The latter have the useful property that they may be numerically simulated, under 
certain conditions, and we provide the Kubo oscillator as a numerical test case. 

In Chapter 3, we present a theory of cascade two-photon emission in an atomic 
ensemble. The radiative atomic dipole-dipole coupling is shown to infiuence the 
emission of the idler photon, resulting in the appearance of superradiant time scales. 
The theory is developed on the basis of Schrodinger probability amplitudes assuming 
single atomic excitations. This approach allows a straightforward treatment of the 
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spectral entanglement properties of the signal-idler photons. 

In Chapter 4, we relax the assumption of single atomic excitations and develop a 
theory based on c-number stochastic partial differential equations, derived using the 
methods reviewed in Chapter 2. Numerical solutions of the equations are used to 
compare with the superradiant timescales derived in the analytical theory. 

In Chapter 5, the analysis of Chapter 3 is used to discuss the behavior of the 
cascade emission on the DLCZ protocol for the quantum repeater. Entanglement 
swapping and quantum teleportation are investigated, and the influence of time- 
frequency entanglement is discussed. 

In Chapter 6, the use of the diamond configuration in frequency up and down con- 
version is analyzed using quantum-Heisenberg Langevin and Maxwell-Bloch equation 
methods. We present results for the optimal conversion efficiency as a function of 
optical thickness of the atomic ensemble. The role of pulse length and quantum 
fluctuations are discussed. 

In Chapter 7, we present some conclusions. 

In Appendixes A-D, we present a great deal of supporting information on the theo- 
retical derivations that are quite lengthy on account of both the multimode treatment 
of the light fields and the complicated atomic level schemes and atomic dipole-dipole 
interactions. 
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CHAPTER II 



REVIEW OF THEORETICAL AND 
NUMERICAL METHODS 

In this Chapter, we review the derivations of quantum-Heisenberg equations and c- 
number Langevin equations for hght-atom interactions. The reason for focusing on 
these methods is, in the first place that they are less familiar than Schrodinger picture 
methods (see Chapter 3 and 5) and that our applications of these methods (Chapter 
4 and 6) involve rather long derivations that may obscure the basic ideas. 

We provide two methods of deriving the c-number Langevin equations and their 
noise correlations. The equations may be found from the quantum Heisenberg- 
Langevin equations using a "quantum-classical" correspondence [26]. Alternatively, 
c-number Langevin equations are deduced by a Schrodinger-picture approach that 
employs characteristic equation and coherent state phase space methods. In the 
final step the Langevin equations are deduced from a Fokker-Planck equation for a 
generalized statistical distribution. Such methods were initially applied in quantum 
laser theory in the 1960's by Haken [27]. The independent derivations will be used 
to check the lengthy derivations involved in the case of cascade emission. 

2. 1 Quantum and C-number Langevin Equations 

Langevin equations were initially derived to describe Brownian motion [28]. A fluc- 
tuating force is used to represent the random impacts of the environment on the 
Brownian particle. A given realization of the Langevin equation involves a trajec- 
tory perturbed by the random force. Ensemble averaging such trajectories provides 
a natural and direct way to investigate the dynamics of the stochastic variables. 
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Figure 2.1: The two-level atomic ensemble interacts with a classical and quantum 
field, (a) An elongated atomic ensemble of length L is excited by a pump field of 
Rabi frequency Via and emits a propagating quantized field denoted by annihilation 
operator . (b) Two-level structure for an atomic ensemble with the ground (|0)) 
and excited (|1)) state. The detuning of the pump field is Ai. 
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In this section, we review quantum and c-number Langevin equation approaches 
for a two-level atomic ensemble interacting with a quantized electromagnetic field. 
As shown in Figure 2.1, the atoms are excited by a pump field of Rabi frequency 
and a propagating quantized field is considered to be emitted along the direction 
z of the ensemble with length L. 

2.1.1 Quantum Heisenberg-Langevin equations 

We consider the Hamiltonian of N two-level (ground and excited states |0), |1)) 
atoms interacting with one pump field and a multimode quantized fields with mode 
annihilation operators ai that satisfy the commutation relation [ci;, a],] = Sw for the 
Ith section along the propagation direction. The propagation length L is discretized 
into 2M + 1 elements [29]. In the electric dipole approximation and rotating wave 
approximation, the interaction is given by —d-E, Appendix B.l. The Hamiltonian H 
includes the free evolution {Hq) of atoms with transition frequency oui, the quantized 
field of central frequency ou, and the dipole interaction (i?/), 

H^Ho + Hi, (2.1) 

M M 

Ho^ Yl ^i^nW + ^ Yl a\{t)ai{t) + hJ2^ii'0'\{t)ai'it) , (2.2) 

l=-M l=-M 1,1' 

M 

Hi^-hY [^Lao\{ty''^"'~'"'^' + 9V2M + lal\{t)di{t)e"''^ + /i.e.] (2.3) 

l=-M 

where aQi{t) = ^^^^ |0)/i(l| . The Rabi frequency fli — dioS{kL)/ {2h) is one- 

rij,=zi 

half the conventional definition. The dipole matrix element dio = (l|(i|0), coupling 
strength g = dioS{k)/h where S{k) = y^huj/2eoV is the electric field per photon, 
and Zp = ^fe' P = ^- The matrix ouu, = Y!^^_m kne'^-^'^-'^'^ / {2M + 1) 

accounts for field propagation by coupling the local mode operators. 

The dynamical equations including dissipation due to spontaneous emission can 
be treated by introducing the reservoir field that interacts with the system [30] . After 
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introducing the coupling to the reservoir, we may write down by inspection the dissi- 
pation terms. We define 701 to be the spontaneous emission rate from |1) — )■ |0). In 
the co-moving frame coordinates z and t = t — z/c, the quantum Heisenberg-Langevin 
equations are 



—5-01 = (iAi - ^)5-oi + i^a{^oo - ^11) + «5'(c^oo - cfii)E^ + Jqi, (2.4) 
Q 

—an = -7oia-ii + iQ^loIi - iQ^^'oi + ig^liE^ - ig*^oiE~ + Tu, (2.5) 

jr-E+ = ^^doi + :FE+, 2.6 
oz c 

where various Langevin noises associated with atomic operators aoi, an and field 
operator E'^ are necessary to preserve equal time commutation relations. The detun- 
ing of the pump field is — uil — ^^i and the slowly-varying operators are defined 
as 5-01 = 4ie-**="^'+*'^"7A/'^, an = aJi/iV^, and E+{z,t) = ^2M + laie'"^-* where we 
let ui — ujl. The time evolution of atomic coherence (aoi) depends on the popu- 
lation difference (aoo — an), and in turn atomic population is infiuenced by atomic 
coherence and the classical and quantized fields. The atomic coherence couples to 
the quantized field along the propagation direction, z. 

The noise operator correlations are related to the dissipation through the fiuctuation- 
dissipation theorem [31, 30]. If we have a quantum Langevin equation for variable 

X 

i:{t) ^ A,{t) + F,{t) (2.7) 

where is so-called the drift term for x, and the corresponding Langevin noise 
operator is F^, the quantum noise correlation functions can be derived from the 
generalized Einstein relation, 

where the bracket denotes the quantum mechanical ensemble average. 
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With the above recipe, we have the non- vanishing normally ordered quantum noise 
correlation function from Eq. (2.6), 

= 701^11- (2.9) 

where (^J^li{t, z)^ii{t' , z')"^ — ^S{t — t')S{z — z') (^D\\,i\J^ and D is also referred to 
as a diffusion matrix element by analogy with classical diffusion processes. 

Even for this relatively simple light-matter interaction, there is no analytical solu- 
tion possible. The c-number Langevin equation approach, below, provides a possible 
way to attack the problem numerically by stochastic simulation and to calculate 
normally-ordered quantities by ensemble averaging, although we will not pursue such 
simulations here. 

2.1.2 C-number Langevin equation 

A c-number Langevin equation approach may be suitable for stochastic simulation 
[32, 28], and utilizes the methods developed by Lax, Louisell, and Haken to describe 
the dynamics of the interaction [27, 26]. Their recipe involves a normal ordering 
procedure and a so-called "quantum-classical correspondence" to derive the c-number 
Langevin equations [26, 31, 33]. The normal ordering chosen is aji, an, aoi, E~ , £"+ 
where the creation operators always appear to the left of the annihilation operators. 
The population operator is put between the atomic coherence operators since it is 
self conjugate. 

The c-number Langevin equations are then derived from Eq. (2.6) by making the 
quantum-classical correspondence that we denote as 

a-Ji CKs, CTii q;4, o-Qi as, E~ E" , E^ E^ . (2.10) 
Similarly for the Langevin noises. 
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where the notation is chosen to facihtate the comparison with an alternative approach 
that we will discuss in the next Section. 

The classical noise correlation functions are also derived from an Einstein relation. 
Consider the c-number Langevin equation for the variables x and y, 

x{t) = A,{t)^F,{t), (2.12) 
Z/(i)=^(i) + F,(i). (2.13) 

From the requirement of equivalent time evolution of normally-ordered operators and 
their c-number counterparts, we have for example 

= ^(xy). (2.14) 

Classical noise correlations can be derived from the quantum ones using 

(F.Fj,) = + (x^) + - (^xAy) - {A,y). (2.15) 

where the quantum and classical noise correlations are formally quite different. For 
non-normally-ordered operators xz, we may use the commutator to substitute that 

{xz) ^{zx) + {[x,z\). (2.16) 

The drift term of the c-number Langevin equations are closely related to the corre- 
sponding term in the quantum Heisenberg-Langevin equations. After the quantum- 
classical correspondence is made, we derive the coupled equations with c-number 

variables {E^ , , a^, a^, a^,) and Langevin noises (^^1,2,3,4,5) that satisfy 
d 

—as = {iAi ^)q;3 + iQa(ao - 0:4) + ig{oio - a^jE'^ + T3, (2-17) 

d 

—a4^^-ry^^j^a4 + iflaC(5-i^*aOi3 + WOi5E'^-ig*a3E + J^4, (2-18) 
^0^5 = {-iAi - -y)ci;5 - i^liao - "4) - i9*{ao - a^jE' + (2.19) 

l-V = + J-.. ^E- = + r,, (2.20) 

oz c oz c 
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The associated non-vanishing diffusion matrix elements, however look quite dif- 
ferent to their quantum counterparts 

-03,3 -i2flaOi3 - i'^'QCtsE^ , 

D4 4 — iQaCKs — i^a^3 + ioi^E'^ — ia^E~ -\- 7010:4. (2-21) 

The diffusion matrix elements are defined as {Ti{t, z)J^j{t' , z')) = §S{t - t')S{z - 
z') (Dij) in the continuous limit. For the more complicated light-matter interactions 
we will encounter in Chapter 4 involving four atomic levels interacting with two prop- 
agating quantized light fields, the diffusion matrix calculation is much more intricate. 
It is therefore important to have an independent check of the c-number equations 
and the associated diffusion matrix. In the following Section we review the Fokker- 
Planck equation approach based on a Schrodinger picture treatment of the quantized 
light-atom interaction. 

2.2 Fokker-Planck Equations and Stochastic Dif- 
ferential Equations 

Here we review the alternative method, due to Haken [27], to derive the c-number 
Langevin equations or equivalently stochastic differential equations via Fokker-Planck 
equations [27, 35, 32, 34]. 

The Fokker-Planck equation is used to describe the fluctuations in Brownian mo- 
tion [28], and its solution for probability distribution f{x,t) of Brownian particles in 
space X and time t is determined by the drift and diffusion properties of the particles. 

2.2.1 Characteristic functions in P-representation 

The Characteristic function x is convenient for the derivation of Fokker-Planck equa- 
tion, and it is the distribution function of the Fokker-Planck equation in Fourier space. 
We follow the same procedure of P-representation laser theory [27] . 
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The relevant operators of our system are atomic coherences (ctq^, ctq^), population 
(crj^) and field operators {a], ai). The normally ordered exponential operator is 
chosen to be 



E{\)^Y[E\\l 

I 

E\X) = e*'^BC^oie*'*'4*iie''^3*oie''^2«(e'^i«i^ (2.22) 

where E[\) the complete exponential operator and is decomposed into products of 
E^{\) for each section I of the propagation direction. We note that the ordering of 
operators is the same as we chose for the quantum-classical correspondence in the 
previous Section. The complex parameters are classical counterparts of operators 
in Fourier space, as will become clear when we derive the Fokker-Planck equation. 
Then characteristic function % can be calculated from a density matrix p, 



X = Tt{E{X)p}, (2.23) 
|=-H|}-E(|) . » = (2.24) 
and time evolution of p is 

where Hq — Ha + H^. Ha is the Hamiltonian for atomic free evolution, Hl is the 
Hamiltonian for the pump field, and the dipole interaction Hamiltonian is Hi. The 
dissipation from spontaneous emission is denoted as sp. 

The contribution from (g^)^^ is calculated up to the second order in Aj. The 
validity of truncation to second order is due to the expansion in the small parameter 
1/Nz. The dissipative contribution, identified by superscript (2), takes the form. 
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7oiTV|£;(A) 

iXs d 



7oi 



(7oipo-oi - -aup - -pau 

— 1X4- 



(2) 



+ 



{ix^y d 



(2.25) 



2 d{iX^) 2 ^(iAs) "d{iXi) 2 ^(^4). 

where we drop the summation over spatial shces Z, which we will retrieve later. Col- 
lecting together all contributions to the characteristic function, we may proceed to 
write down a Fokker-Planck equation that leads to the c-number Langevin equation. 

2.2.2 A Complimentary Derivation of C-number Langevin Equations 

The time derivative of the distribution function / is found from the Fourier trans- 
form of the characteristic function ^ = j^^^ / ••• / e~'^°''^^dXi...dXn- Separating the 
different contributions we may write 



dl 
dt 



1,1' 



(2.26) 



The details of the C operators can be found in Appendix B. Here we show Cj as 
an example, 



an) - 



_d_ 
dai 



;-2ai + 7V,) + e ^(4) 



+ igV2M + le'''^' 



+ ^g*^/2M + le-'''^{4) (^4 - ^) + {C'T , 



(2.27) 



where C" is the correspondence that -H- a^^ a\ -H- a^^ a[ •<->■ 0:2, and * denotes 
complex conjugation. The results is a Fokker-Planck equation of the form 



df 



d 



d 



dt da"^"-^ 813^^^ ^ 2\dadl5 ' d(3da 



(2.28) 
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where A^^is and D^^ are drift and diffusion terms. The corresponding c-number 
Langevin equations may be derived rigorously when D is positive definite, and take 
the form 



^^A^ + T^,^=A^ + r^ (2.29) 

with a classical noise correlation (rQ(i)r^(i')) — S{t — t')Dap. Higher order deriva- 

a_ 

tivcs (third order and higher, from the Taylor expansions of e ) are ignored as 
they involve the small parameter l/N^. The corresponding c-number Langevin, or 
stochastic differential, equations are 



+ ig\/2M + le'^'' (q;|, - a\)a[ + T^, (2.30) 

+ igV2M + le'^''a[a\ - igW2M + le-'''''a[a[ + r[, (2.31) 

a{ = -iuja\ - i^0Jwa{ + V2M + le'^^^'^ + (2-32) 
I' 

We can retrieve the continuous limit with the slowly varying variables, a^{z,t) = 
Q,^g-ifc„^,+ia;aY^^^ Q;4(^,t) = a\/N^, E+{z,t) = V2M + la^e^'^"*, and note that 

—i Uiiia^ = —c-^a\ and ctg = A^^ — ol'/^. Define also the slowly-varying 
Langevin noises, 

J'i(^,i) = \/2M + le''"*r'i. (2.33) 
Finally, in the co-moving frame coordinates z and r = t — z/c, the c-number 
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Langevin equation becomes 



d 



(iAi z-)a;3 + i^a{oiQ - a.^) + ig{a.Q - ai)E^ + J's 



(2.34) 




-7oi«4 + i^a^^ - i^lois + iga^E^ - ig*azE + J4, 



(-iAi z-)q;5 - ifia(Q;o - ^a) - ig*{ao - a4)E + T5, 



iNg* d '^9 

c oz c 



3, 



(2.36) 



(2.35) 



(2.37) 



where Ai — oua — uji. The non- vanishing diffusion coefficients extracted from the 
Fokker- Planck equation are 

-03,3 = -i2D,aa3 - i2a^Ep, 04^4 = if^aSs - i^lois + ia^E^ - ia^E^ + 7oi54. (2.38) 

Comparing with the results in the previous Section and Eq. (2.21), we find com- 
plete agreement. As the c-number Langevin equations are derived from a Fokker- 
Planck equation, they should be interpreted as Ito-type stochastic differential equa- 
tions (SDE), and this is important in the numerical solution method [32]. In nu- 
merical simulation it is common to first transform from the Ito equation to its corre- 
sponding Stratonovich form. 



We present an example of the Kubo oscillator to illustrate numerical simulation of 
a multiplicative noise stochastic differential equation. A Kubo oscillator provides 
a good test case in the numerical solution of stochastic differential equations. The 
Langevin equation of the dimensionlcss Kubo oscillator with amplitude z{t) is given 
by the Stratonovich equation. 



where ^{t) is a delta-correlated real Gaussian distributed noise with zero mean, 
{C{t)) = , and {i{t)^{t')) = 6{t — t'). The bracket denotes an ensemble average. The 



2.3 Kubo Oscillator 



d 

dt 



z{t) = z^{t)z{t) 



(2.39) 
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1 2 3 4 5 

t (arb. units) 



Figure 2.2: Kubo oscillator simulation. The time evolution of Re(z(t)) (dashed- 
red) averaged from an ensemble of 1024 simulations. z{0) = 1. We compare with the 
exact solution, z(0)e~*''^ (solid-black), and find good agreement. A demonstration 
of one stochastic realization (dashed-circle blue) shows large fluctuation around the 
averaged and exact results. Note that the imaginary part of the solution is almost 
vanishing as it should be, and is not shown here. 

exact analytical solution for the first moment is {z{t)) = {z{t = 0))e~*/^. To numer- 
ically simulate the Stratonovich equation (2.39), we use the following discretization 
in time [36, 37] 

Z{tm) = Zn-1 + i^{tn-l)z{tm)^ (2.40) 

where z{tm) is evaluated at the midpoint, tm = {tn + in-i)/2 and At = t„ — is 
the time step. In this specific case where the noise is linear in z(t), we may solve Eq. 
(2.40) to give z{tm) = Zn-i/ (1 - i^(t„-i)At/2). Setting z{tn+i) = 2z{tm) - z{tn), we 
use z{tn+i) for the next time step of the integration. 
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The Langevin noise is numerically simulated as ^{t) — randn(i)/-\/Af , where 
randn(t) is a random number generated from a Gaussian distribution with zero mean 
and unit variance. In Figure 2.2, we compare the analytical and numerical results for 
the Kubo oscillator. The initial condition is set as z{0) — 1, and we use 1024 realiza- 
tions for the converged numerical result {z{t)) with At = 0.01. The numerical result 
is in good agreement with the exact solution, 2;(0)e~*/^. The temporal evolution of 
one typical realization of the stochastic process fluctuates significantly, as shown. 

We will use the approach demonstrated here to simulate the more complicated 
c-number Langevin noises in our investigation of cascade emission from an atomic 
ensemble in Chapter 4 (see also Appendix B). 
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CHAPTER III 



SUPERRADIANT EMISSION FROM A 
CASCADE ATOMIC ENSEMBLE: 
ANALYTICAL METHOD 

In this Chapter, we use Schrodinger's equation to investigate cascade emission from a 
four-level atomic ensemble. Quantum communication has opened up the possibility 
to transmit quantum information over long distance. Due to the transmission loss 
in long distance fiber-based quantum communication, telecommunication (telecom) 
wavelength light is important to maximize the transmission efficiency. The alkali 
atomic cascade transition shown in Figure 3.3 is able to generate telecom wavelength 
light, the signal, from the upper transition and a near-infrared field, the idler, from the 
lower one. The telecom light can travel through the fiber with minimal loss, while the 
near-infrared field is suitable for storage and retrieval in an atomic quantum memory 
element. Their use in a quantum information system requires quantum correlations 
between stored excitations and the telecom field. 

We develop a quantum theory to characterize the properties of the correlated sig- 
nal and idler photons and study how the laser excitation pulse modifies their spectral 
profile. The wave packets of this entangled source are found, and Schmidt decom- 
position provides the basis for engineering a pure photon source that is crucial in 
quantum information processing. 



21 



3. 1 Introduction 

The spontaneous emission from an optically dense atomic ensemble is a many-body 
problem due to the radiative coupling between atoms. This coupling is responsible 
for the phenomenon of superradiancc firstly discussed by Dicke [23] in 1954. 

Since then, this collective emission has been extensively studied in two atom sys- 
tems indicating a dipole-dipole interaction [20, 21], in the totally inverted N atom 
systems [38, 39], and in the extended atomic ensemble [22]. The emission inten- 
sity has been investigated using the master equation approach [40, 41, 42] and with 
Maxwell-Bloch equations [43, 44]. A useful summary and review of superradiancc 
can be found in the reference [45, 46]. Recent approaches to superradiancc include 
the quantum trajectory method [47, 48] and the quantum correction method [49] . 

In the limit of single atomic excitation, superradiant emission characteristics have 
been discussed in the reference [50] and [51]. For a singly excited system, the basis 
set reduces to N rather than 2^ states. Radiative phenomena have been investigated 
using dynamical methods [52, 53, 54] and by the numerical solution of an eigenvalue 
problem [55, 56, 57, 58]. A collective frequency shift [59, 60] can be significant at 
a high atomic density [61] and has been observed recently in an experiment where 
atoms are resonant with a planar cavity [62] . 

3.2 The example of two- state atoms interacting 
with a pump field 

The atomic dynamics of two-state atoms interacting with a pump field generally 
requires a basis of 2^ orthogonal states. In this Section we investigate multiple ex- 
citations by a laser by solving numerically the master equation for few atom systems 
(A" = 2, 3, 4), using the quantum optics toolbox [63]. The complete orthogonal states 
may be chosen as 1 symmetric state and {C^ — 1) non-symmetric states for any excita- 
tion number n where is the combination coefficient. It is natural to construct the 
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complete orthogonal states using this decomposition because the interaction Hamil- 
tonian of the pump field, Hi = [-h^ 1 1) (0 1 e^^'^^*^ + c.c] - ^Ai |1)(1|, has the 
same form for each atom. 

For the example of two two-state atoms, there are 4 orthogonal basis states: the 
ground state |00), the symmetric state of a sing le excitation (e^*^-*'i |10) + e*^-^"'2 |01))/v/2, 
the associated anti-symmetric state (e'^^'^^llO) — e'^'^'^^\01))/^/2, and the state of two 
excitations e^'^'^^'i+^^^lll). Note that the spatial phase factor for different atomic 

— * 

position r is included due to the pump field of the wavevector k that is directed along 
the z axis. If more atoms are involved, the complete states of multiple excitations 
can be derived by extending the results of reference [52], and here we list the states 
of four atoms {N — A), 



n = 0, |0i) = |0, 0,0,0), 

I = 1,2, ...AT - 1; 



n = 1 < 



l)iv|0) 
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n = 2 < 



^*''"-^^'"-|l)^-i|l)A.|0)A^iv-i,iv, m = 1,2,...7V2 - 1, 



/N2 



where N2 = N{N -l)/2; 



N 



(3.1) 



where (j, /) in the subscript of the Kronecker delta function of two excitation states 
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is defined so that (1, 2) = 1, (1, 3) = 2, (1, 4) = 3, (2, 3) = 4, and (2, 4) = 5. Note that 
the n = 0, \4>i), and n = N, |02^)) states are symmetric. For n ^ 0,N excitations, 
the states are constructed from one symmetric state and non-symmetric states. 

To investigate the probabihty of multiple atomic excitations in conditions of weak 
off-resonant excitation, we choose a configuration of four atoms that sit on the vertices 
of a square with side dg. The atomic density matrix includes a laser excitation term 
in addition to the one and two-atom dissipation terms; these arise from spontaneous 
emission and radiative coupling due to dipole-dipole interaction [21]; see Eq. (A. 15) 
in Appendix A. We numerically solve for the time evolution of the density matrix. 
The result of steady state single- and double-excitation populations are shown in 
Figure 3.1 as a function of dg. We have assumed a continuous laser field with peak 
Rabi frequency fla = 0.2j and detuning Ai = 57, where 7 is the single-atom sponta- 
neous decay rate for the excited state. The populations of the symmetric states are 

=Tr(p|02)(02|) for a single excitation and P| =Tr(/}|</)6)(06|) for double excitations 
where p is the density operator of the atomic system. The total populations of the 
non-symmetric excitation states are =Tr(p(^^^2 , for a single exci- 

tation, and P"^ =Tr(p (XlxLy \4>x){4>x\)) , for double excitations, respectively. The 
probabilities of three and four excitations are negligible under the weak excitation 
conditions we consider. 

As dg approaches and exceeds A (the transition wavelength), the populations tend 
to the independent atom limit when dipole-dipole coupling is omitted. In this limit, 
the probability of exciting any non-symmetric states goes to zero. The single and 
double excitation probabilities, P^ and P|, are normalized to their independent atom 
values, P^°^ = PeP^Cf and Pf ^ = PePgCl where Pe = nl/{4Al + 7^) [64], Pg = 
1 — Pe, and C" is the combination coefficient. For dg <^ X, the populations of the non- 
symmetric states are comparable to the symmetric ones, indicating the importance of 
dipole-dipole interactions. We see no evidence of a dipole blockade effect in this limit 
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Figure 3.1: Single- and double-excitation populations as a function of distance 
dg. (a) The populations of the symmetric state for a single excitation PI (dashed- 
red) and the sum of non-symmetric single-excitation states Pi'^ (dashed-dotted black), 
(b) The populations of the symmetric state for double excitations P| (dashed-blue) 
and the sum of non-symmetric double-excitation states P^^ (dashed-dotted black). 
P^' and P|' are normalized respectively by the solutions of non-interacting atoms 
Pf ^(solid-red) and Pf ^ (solid-blue). P}°^ = 1.58 x lO'^ and Pf ^ = 9.4 x 10"^ are 
the single- and double-excitation probabilities for independent atoms. 
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for four atoms, but we have observed it in the case of two atoms. Dipole blockade 
refers to the predominance of single excitations as dipole shifts detune double and 
higher excitation states. 

In Figure 3.2, we show the time evolution of -Pf(t) and -P|(i) for dg — 3A (this 
corresponds to an atomic density 8 x 10^° cm~^). The period of the Rabi oscillation is 
determined by 27r/Ai, and the asymptotic steady state value for P| is about 1.6x10"^. 
This coincides with the approximate result |-\/]VQa/(2Ai)p that is found when we 
truncate the basis to the ground state and the orthogonal states of a single atomic 
excitation. 

We also numerically solve a line of atoms (A?" = 2, 3, or 4) with an equal separation 
from dg — Ito 5A, and the results of steady state populations indicate the condition for 
truncation of the basis set at a single atomic excitation is valid when | Ai| \/NQ,a/2. 
If the condition of a single atomic excitation Ai » y/Nfla/'^ is relaxed, we will also 
have dynamical couplings between symmetric and non-symmetric states (at least for 
ds < 3A). It is the dipole-dipole interaction that couples the non-symmetric and 
symmetric states in the presence of the pump laser. 

3.3 Theory of Cascade Emission 

We consider N cold atoms that are initially prepared in the ground state interacting 
with four independent electromagnetic fields. As shown in Figure 3.3, two driving 
lasers (of Rabi frequencies fla and fib) excite a ladder configuration |0) |1) |2). 
Two quantum fields, signal ds and idler Oj, are generated spontaneously. The atoms 
adiabatically follow the two excitation pulses and decay through the cascade emission 
of signal and idler photons. Based on the discussion in the previous Section, we permit 
only single atomic excitations under the condition of large detuning, Ai » -\/]VQa/2. 
The Hamiltonian and the coupled equations of the atomic dynamics are detailed in 
Appendix A. 
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Figure 3.2: Time evolution of populations for symmetric states Pf and P|. The 
population of the symmetric state for a single excitation is Pf (dashed- red), and that 
for the symmetric state for double excitations is P| (dashed-dotted blue). The pump 
condition is the same as in Figure 3.1 for ci^ = 3A. 
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Figure 3.3: Four-level atomic ensemble interacting with two driving lasers (solid) 
with Rabi frequencies fla and fib. Signal and idler fields are labelled by and dj, 
respectively and Ai and A2 are one and two-photon laser detunings. 

To correctly describe the frequency shifts arising from dipole-dipole interactions, 
we do not make the rotating wave approximation on the electric dipole interaction 
Hamiltonian. The frequency shift has contributions from the single atom Lamb shift 
and a collective frequency shift. The Lamb shift is assumed to be renormalized into 
the single atom transition frequency distinguishing it from the collective shift due to 
the atom-atom interaction. 

3.3.1 Probability amplitudes for signal and signal-idler emissions 

Writing the state- vector {ipit)) in a basis restricted to single atomic excitations, and 
single pairs of signal and idler photons, we can introduce the probability amplitudes, 

N 

Cs,kXt) = ^e-*-''(3„l,,,,J^(t)) (3.2) 

and 
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^mW = (0,1/c3,a.,1..aJ^W) (3-3) 
defined in Appendix A. Note that Cs^kXt) is an ampfitude for a phased excitation of 
the ensemble of atoms subsequent to signal photon emission. 

After adiabatically eliminating the laser excited levels in the equations of motion, 
we are able to simplify and derive the amplitude Cs,ki and the signal-idler (two- 
photon) state amplitude Ds,i as shown in Appendix A, 

Jo 



Jo Jo 



Jo Jo 

gi(a;i-a;3)t'gi(a;.-a,23-A2)t"^(^//)_ (35) 

The factor Q^^k-Tf, reflects phase-matching of the interaction under conditions 
of four-wave mixing when the wavevector mismatch Ak = ka + ki, — ks — ki ^ 0. The 
radiative coupling between atoms results in the appearance of the superradiant decay 
constant 



r^^{Nfi + i)rs (3.6) 

where is the natural decay rate of the |3) — )■ |0) transition, and /i is a geometrical 
constant depending on the shape of the atomic ensemble. An expression for the 
collective frequency shift Soui is given in the Appendix A. As shown in Figure 3.4, 
we numerically calculate the geometrical factor /2, Eq.(A.12), to demonstrate how the 
decay factor Np, + 1 depends on the height and radius of a cylindrical ensemble. The 
arrows in the figure point out the contour lines (yellow and green) of Np, + 1 A and 
6 which are comparable to the operating conditions of the experiment [16]. 
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Figure 3.4: The superradiance decay factor (/i = fi) for a cylindrical ensemble 

of length h and radius a in unit of transition wavelength A. The atomic density is 
8 X 10^° cm~^ and A = 795 nm corresponding to the Dl line of ^^Rb. See the text 
for the explanation of the arrows. 
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In the above expressions for the probabihty amphtudes, b{t) — "^^^^^ is propor- 
tional to the product of the Rabi frequencies. We use normahzed Gaussian pulses 
as an example where flait) — -^^^a'^~^^ ^b{t) = -^^^b^'^^^^^ , so that the two 
pulses are overlapped with the same pulse width. Cla,,b is the pulse area, and let 
Acus = cug — UJ23 — A2 — SijJi, Aui = uji — u!3 + SuJi. We have the probability amplitude 
for signal photon emission and atoms in a phased state, 



Cs,ki{t,Au}s) 



4A1A2 



E 



4A1A2 ^ TTT^ 2 V 2 



pJV T,JV 



and the two-photon probability amplitude is 

Dsi{t, Aujs, Auji) 



2v/2r 

+ ,-(A...A.,)V^,f . (^1 + erf( 4^-^(^^^+^"-)"^ )) }, (3.8) 
where erf is the error function 

erf(a;) = ^ / e~^'dt. (3.9) 

VT^ Jo 

Asymptotically Dsi approaches the value. 



DJAco,, Acui) = \ . ^ '1^^J±^ , (3.10) 
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indicating a spectral width for idler photon in a Lorentzian distribution mod- 

ulating a Gaussian profile with a spectral width 2\/2/r for signal and idler. Energy 
conservation of signal and idler photons with driving fields at their central frequen- 
cies corresponds to ujs -\- ijJi — ijJa i^hi which makes AO's -|- Awj = 0; the collective 
frequency shifts cancel. 

3.4 A Correlated Two-photon State 

Using the asymptotic form of the two-photon state given in Eq. (3.10), the second- 
order correlation function G\ / is calculated as [30] 



= m ^ ^)\E:{T,M)k{r2MEt{r2M)Et{riM)m ^ 00)) = |$,,|^ 



(3.11) 

= {0\E+{f,,h)E+{n,t,)m ^ 00)) (3.12) 



hs 5 A 



EHr2,h) = Y.\J ^^^^k.^A''-'^''' (3.14) 

/cjjA 

where \ip{t 00)) denotes the state vector in the long time limit that involves the 
ground state and two-photon state vectors. Free electromagnetic fields, signal and 
idler photons, at space {f\,r2) and time (^1,^2) are and E^ where (-I-) denotes 
their positive frequency part. For second order correlation function, only Dgi, derived 
in the previous Section contributes to it, then we have. 
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2 « 

ci0s[4 - ks{ks ■ ds)] 



4A1A2 4e§c6(27r)6 

gi(A.,-5..)(ti-^) y dA^^e-'^'^^(*i-^)e-^^»'^'/« (3.15) 

where we have used the change of variables in the first step, replaced (jJs = UJ23 + A2 + 
Aujs + 5uji, and changed the variable AO's — >■ Aoug — Aoui. Solid angle integration is 
denoted as dQs,i for signal (idler) photon. The divergent part of (which varies 
relatively slowly) has been moved out from the integral of dAcUs, and we replace cUg 
with the signal transition frequency 0023- We then have 



iAk-fu, I 7* 1 1 7* I .3 3 



4A1A2 4egc6(27r)6 
r J (L 



dflsdfli[ds - ks{ks ■ ds)][di - ki{ki ■ di)] 



— iAui — iSuji) 

(3.16) 

where we replace uji — 003 + Aoui — Soui and change the variable Aoui — > Auji + Soui. 
The divergent part of cuf is again moved out from the integral of dAoui and replace uji 
with the signal transition frequency (jj^. Finally we have 
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= I^^^^^^'^"^'''' I - ^^^^^ • - ^^^^^ • ^'^^ 

fj, 

e(t,-^-ti + ^) (3.17) 

c c 

where the complex integral with the pole at Acjj = in the lower half 

plane leads to a step function © that shows the causal connection between signal and 
idler emission. The emission time for the signal field {ti — ^^^) is within the pulse 
envelope of width r, and the idler photon decays with a superradiant constant 
Note that the collective frequency shift 5u!i appears in the signal {0023 + A2 + Suji) and 
idler (us — Sui) frequency consistent with energy conservation. Let Atg = h — 
and Ati = ^2 - we then have 



= ^ip^^ S^^^''"'' / ^^^^^^^^"^ - ^^^^^ • ^^^^^^^ - • ^^^^ 

If we let At = Ati — and choose Atg = as the origin in time (idler gating 
time) , then we have the second-order correlation function 

G^JiAt) = |$.,i(At)|2 oc e-^3"'At ^here At > 0. (3.19) 

It resembles the result for the second-order correlation function in the case of single 
atom, whereas here we have an enhanced decay rate due to the atomic dipole-dipole 
interaction. 
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In Figure 3.5, we plot out the absolute value of spectrum Dgi{AuJs, Awj) and the 
second-order correlation function G). l{Ats, Ati). In (c), we show for T^Ati = 0.2. 
The width of 1/Ati = bF^ corresponds to = {Njl + 1)F3 = bF^. 

3.5 Schmidt Decomposition 

Correlated photon pairs may be generated by parametric down conversion (PDC) 
[65, 66, 67]. The degree of entanglement can be quantified by Schmidt mode decom- 
position [68, 69], allowing the influence of group- velocity matching [70] to be assessed. 
A pure single photon source is a basis element for quantum computation by linear 
optics (LOQC) [71], and it can be conditionally generated by measurement [72]. A 
similar approach can be applied to the study of the transverse degrees of freedom 
in type-II PDC [73] and PDC in a distributed microcavity [74]. In photonic-crystal 
fiber (PCF), a factorizable photon pair can be generated by spectral engineering [75]. 
The spectral effect has been discussed in relation to a quantum teleportation protocol 
[76] as a first step toward quantum communication. 

We would like to perform an analysis of entanglement properties of our cascade 
emission source. In addition to polarization entanglement, a characterization of 
frequency space entanglement is required to clarify its suitability in, for example, the 
DLCZ protocol [6]. 

In the long time limit, the state function is given by, Eq. (3.10), 

\^) = |0, vac) +5] iP^^IO, 1,-^^,^,1,-^,; (3.20) 

s.i 

where s — {ks,Xs), i — {ki,Xi), and |0,vac) is the joint atomic ground and photon 
vacuum state. 

The spatial correlation of two-photon state in FWM condition can be eliminated 
by pinholes or by coupling to single mode fiber so we consider only the continuous 
frequency space. For some specific polarizations and Aj, we have the state vector 
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Figure 3.5: (a) Absolute value of the spectrum for two-photon state probability 
amplitude Dg i and (b) the second-order correlation function G), ({Atg, Ati). (c) A 

normalized G].l{Ats = 0, Atj) with T^r = 0.2. The exponential decay corresponds 
to the superradiant decay factor Nfi + 1 = 5. 
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I*) = J /K,c^i)ai^(t^,)ai^(c^,)|0)(iw,(ia;i, (3.21) 

where 

f(Us, Ui) = — , . (3.22) 



— iAcoi 



The quantification of entanglement can be determined in the Schmidt basis where 
the state vector is expressed as 



l*)=Ev^^n^n|0)' (3.23) 

n 

bl = J ^n{cOs)alMdWs, (3.24) 

ci = J (f)n{i^i)a{.{u;i)duji, (3.25) 

where c]^ are effective creation operators. Eigenvalues A„, and eigenfunctions tpn 
and 0„, are the solutions of the eigenvalue equations, 



Ki{uj,Uj')'ll!n{^^')duj' = Xnlpni^^), (3.26) 

K2{uj,uj')(t)n{uj')(ko' = A„0„(a;), (3.27) 



/ 



where Ki{(jj,u') = J f{u,Ui)f*{u',Ui)du!i and K2{u,u') = f f{u2,u;)f*{u2,uj')du2 
are the kernels for the one-photon spectral correlations [68, 69]. Orthogonality of 
eigenfunctions is J ipi{uj)ipj{uj)duj — 5ij, J (f)i{uj)(l)j{uj)duj — 5ij, and the normalization 
of quantum state requires A„ = 1. 

In the Schmidt basis, the von Neumann entropy may be written 

oo 

S^-Y, AnlnA„. (3.28) 

n=l 
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If there is only one non-zero Schmidt number Ai = 1, the entropy is zero, which 
means no entanglement and a factorizable state. For more than one non-zero Schmidt 
number, the entropy is larger than zero and bipartite entanglement is present. 

The kernel in Eq. (3.22) has all the frequency entanglement information, entan- 
glement means /{cug, uji) cannot be factorized in the form g{uJs)h{(xJi), a multiplication 
of two separate spectral functions. By inspection the Gaussian profile of signal and 
idler emission is a source of correlation. The joint spectrum Aoug + A.u!i is confined 
within the width of order of 1/r. The Lorentzian factor associated with the idler 
emission has a width governed by the superradiant decay rate. 

In Figure 3.6, we show the Schmidt decomposition of the spectrum. We use a 
moderate superradiant decay constant Np, -|- 1 = 5, comparable to the reference [16], 
and a nanosecond pulse duration r = 0.25 (Ri 26/4 ns), and r3/27r = 6 MHz. Due to 
slow convergence associated with the Lorentzian profile, we use a frequency range up 
to ±1200 (in unit of F^) with 2000 x 2000 grid. The numerical error in the eigenvalue 
calculation is estimated to be about 1% error. In this case, the largest Schmidt 
number is 0.8 and corresponding signal mode function has a FWHM Gaussian profile 
4-^2 ln(2)/r ~ lOFs. The idler mode function 0i reflects the Lorentzian proflle in 
the spectrum at the signal peak frequency {Aoug — 0), 

/(A..=0,A..) = ^^^^;^55^ (3.29) 

where a relatively broad Gaussian distribution is overlapped with a narrow spread of 
superradiant decay rate [FWHM > {Njl+ l)r3/2]. 

Figure 3.7 shows that the cascade emission source is more entangled if the super- 
radiant decay constant, or the pulse duration increases. We note that the Gaussian 
profile aligns the spectrum along the axis AoUg = — AcUj and the spectral width for sig- 
nal photon at the center of the idler frequency distribution {Acui — 0) is determined by 
pulse duration r. For a shorter pulse > {Nj2 + l)rs/2, the joint Gaussian profile 
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I 3 I 3 

Figure 3.6: Schmidt mode analysis with pulse width r = 0.25 and superradiance 
decay factor Nfi + 1 = 5. (a) Schmidt number and (b) signal mode functions: Re['?/'i] 
(solid- red) and Re[V^2] (solid-blue). Imaginary parts are not shown, then are zero, (c) 
Real (solid) and imaginery (dotted) parts of first (red) and second (blue) idler mode 
functions, 0i and 02- (d) The absolute spectrum \f{Aus,AcOi)\. 
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has a larger width, and the spectrum is cut off by the Lorentzian idler distribution. 
A larger width leads to a less entangled source and distributes the spectral weight 
mainly along the crossed axes AoUg — and Acui — 0. A narrow Lorentzian profile 
cuts off the entanglement source term ^-i'^'^^+^^i)'''^'^/^ tilting the spectrum along the 
line Aujs + Acui = 0. In the opposite limit, < {Np, + l)rs/2, the spectrum is highly 
entangled corresponding to tight alignment along the axis Acug — —Acui (Figure 3.7 
(c)). 

Note that the short pulse duration (r > 0.25 (6.5 ns)) should not violate the 
assumption of adiabaticity r > 1/Ai or I/A2. 

The Schmidt analysis and calculation of von Neumann entropy shows that signal- 
idler fields are more entangled if the ensemble is more optically dense, corresponding 
to stronger superradiance. For the DLCZ protocol, we wish to avoid frequency 
entanglement. The superradiance may be reduced with smaller atomic densities but 
good qubit storage and retrieval efficiency require a moderate optical thickness [16] . 
A better approach involves using short pulse excitation > {Np, + IjPs. We will 
investigate the spectral properties in more details for the DLCZ scheme in Chapter 
5. 
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Figure 3.7: Absolute spectrum of two-photon state and the eigenvalues of Schmidt 
decomposition. Np, + 1 = 5 for both (a) r = 0.25 (b) r = 0.5. Nji + 1 = 10 for (c) 
r = 0.25. The von Neumann entropy (S) is indicated in the plots. 
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CHAPTER IV 



SUPERRADIANT EMISSION FROM A 
CASCADE ATOMIC ENSEMBLE: NUMERICAL 

APPROACH 

In this Chapter, we investigate the cascade emission (signal and idler) from an atomic 
ensemble using a numerical approach. In Chapter 3, we studied the correlated emis- 
sion using Schrodinger's equation assuming single atomic excitations. To relax the 
assumption of single atomic excitations, we derive a set of c-number stochastic differ- 
ential equations derived using the quantum statistical methods reviewed in Chapter 
2. We solve numerically for the dynamics of the atoms and counter-propagating 
signal and idler fields. The signal and idler field intensities are calculated, and the 
signal-idler correlation function is studied for different optical depths of the atomic 
ensemble, and compared with the analytical results of Chapter 3. 

4 . 1 Introduction 

To account for multiple atomic excitations in the signal-idler emission from a cascade 
atomic ensemble, the Schrodinger's equation approach becomes cumbersome. An 
alternative theory based on c-number Langevin equations as discussed in Chapter 
2, is suitable for solution by stochastic simulations. An essential element in the 
stochastic simulations is a proper characterization of the Langevin noises. These 
represent the quantum fluctuations responsible for the initiation of the spontaneous 
emission from the inverted [44, 77, 78, 79], or pumped atomic system [81, 80] as in 
our case. 
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The positive-P phase space method [32, 29, 83, 85, 84, 82, 86] is employed to derive 
the Fokker-Planck equations that lead directly to the c-number Langevin equations. 
The classical noise correlation functions, equivalently diffusion coefficients, are al- 
ternatively confirmed by use of the Einstein relations reviewed in Chapter 2. The 
c-number Langevin equations correspond to Ito-type stochastic differential equations 
that may be simulated numerically. The noise correlations can be represented either 
by using a square [87] or a non-square "square root" diffusion matrix [84]. The 
approach enables us to calculate normally-ordered quantities, signal-idler field inten- 
sities, and the second-order correlation function. The numerical approach involves a 
semi-implicit difference algorithm and shooting method [88] to integrate the stochastic 
"Maxwell-Bloch" equations. 

Recently a new positive-P phase space method involving a stochastic gauge func- 
tion [89] has been developed. This approach has an improved treatment of sam- 
pling errors and boundary errors in the treatment of quantum anharmonic oscillators 
[90, 91]. It has also been applied to a many-body system of bosons [92] and fermions 
[93]. In this Chapter, we follow the traditional positive-P representation method 
[94]. 

4'2 Theory of Cascade emission 

The complete derivation of the c-number Langevin equations for cascade emission 
from the four-level atomic ensemble is described in detail in Appendix B. After 
setting up the Hamiltonian, we follow the standard procedure to construct the char- 
acteristic functions [27] in Appendix B.2 using the positive-P representation [32]. In 
Appendix B.3.1, the Fokker-Planck equation is found by directly Fourier transforming 
the characteristic functions, and making a l/N^ expansion. 

Finally the Ito stochastic differential equations are written down from inspection 
of the first-order derivative (drift term) and second-order derivative (diffusion term) in 
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the Fokker-Planck equation. The equations are then written in dimensionless form by 
introducing the Arecchi-Courtens cooperation units [115] in Appendix B.3.2. From 
Eq. (B.73) and the field equations that follow, these c-number Langevin equations in 
a co-moving frame are, 



— TToi = (iAi - ^)7roi + if^a(7roo - ttu) + i9-Itio2 - iA^^t + -^01 (I), 
d 

— 7ro2 = (iA2 - y )7ro2 - i^a^^U + i^bT^Ol + ITTosEfe''"^''^ - i7r32-B+ + J'o2, 

d 

^TTii = -7oi7rii + 7127122 + i^a^^ll - i^l^^Ol - i^b^^u + «^^67I"12 + ^11, 

d 

_7r22 = -727r22 + inbTil^ - iniTri2 + inl^E^e-'^''' - i^iE; e'^''' + T22, 
^7133 = -7o37r33 + 7327r22 - inl^E+ e''^'"' + in32E;e'^''' + ml^E+ - m^^Er 
+ -^33, 

|-7ri3 = -(^Al + I5i±^)7ri3 - ^^^>03 - + iT^,2E: e''^^^ + ii^l^E^ 

d 

^TTos = — ^7ro3 - if^aTTis + i'Ko2E~ e"^'''' + iinoo - 7733)^;+ + J"o3, 

d 703 + 72 , -o t V ^7-+ -iAkz ■ 77.- , X- 

— 7^32 = ^A2 ^ 7r32 + zi267r|3 -z(7r22 -7r33)£;7e - ^7^o2-£/■ +^32, 

(4.1) 

where (I) stands for Ito type SDE. nij is the stochastic variable that corresponds to 
the atomic populations of state \i) when i = j and to atomic coherence when i 7^ j, 
and J^ij are c-number Langevin noises. The remaining equations of motion, which 
close the set, can be found by replacing the above classical variables, Tr*^ — > Trjj., 
i^k)* ^ ^jk, «^)* ^ i^s,iT ^ <i ' and J-;, ^ Jj,. Note that the atomic 
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populations satisfy Tr*j — TTjj. The superscripts, dagger (f) for atomic variables and 
(— ) for field variables, denote the independent variables, which is a feature of the 
positive-P representation: there are double dimension spaces for each variable. These 
variables are complex conjugate to each other when ensemble averages are taken, for 
example {ttju) = (j^]kj i-^ti) = (-^m)* ■ doubled spaces allow the variables 
to explore trajectories outside the classical phase space. 

Before going further to discuss the numerical solution of the SDE, we point out 
that the diffusion matrix elements have been computed using Fokker-Planck equa- 
tions and by the Einstein relations described in Appendix B.3.3. This provides the 
important check on the lengthy derivations of the diffusion matrix elements we need 
for the simulations. 

The next step is to find expressions for the Langevin noises, and the details are 
given in Appendix B.3.4 in terms of a non-square matrix B [35, 84]. The matrix B is 
used to construct the symmetric diffusion matrix D{a) = B{a)B'^{a) for a Ito SDE, 

dxi = A{t, -^t)dt + Bijit. '^t)dWi{t) (I) (4.2) 

where $,idt = dWi{t) (Wiener process) and {^i{t)^j{t')) = Sij5{t - t'). Note that 
B — > BS, where S is an orthogonal matrix (SS'^ — /), leaves D unchanged, so B is 
not unique. We could also construct a square matrix representation B [32, 28, 87]. 
This involves a procedure of matrix decomposition into a product of lower and upper 
triangular matrix factors. A Cholesky decomposition can be used to determine the B 
matrix elements successively row by row. The downside of this procedure is that the 
B matrix elements must be differentiated in converting the Ito SDE to its equivalent 
Stratonovich form for numerical solution. 

The Stratonovich SDE is necessary for the stability and the convergence of semi- 
implicit methods. Because of the analytic difficulties in transforming to the Stratonovich 
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form, we use instead the non-square form of B [84] that is shown exphcitly in Ap- 
pendix B.3.4. 

In this case a typical B matrix element is a sum of terms, each one of which is a 
product of the square root of a diffusion matrix element with a unit strength real (if 
the diffusion matrix element is diagonal) or complex (if the diffusion matrix element 
is off-diagonal) Gaussian unit white noise. It is straightforward to check that a B 
matrix constructed in this way reproduces the required diffusion matrix D — BB^. 

As pointed out in the reference [86], the transverse dipole-dipole interaction can 
be neglected and nonparaxial spontaneous decay rate can be accounted for by a single 
atom decay rate one if the atomic density is not too high. We are interested here 
in conditions where the ensemble length L is significant and propagation effects are 
non-negligible, and the average distance between atoms d — ^/V/N is larger than 
the transition wavelength A. The length scales satisfy A < d <^ L, and we consider a 
pencil-like cylindrical atomic ensemble. The paraxial or one-dimensional assumption 
for field propagation is then valid, and the transverse dipole-dipole interaction is not 
important for the atomic density we focus here. 

4-3 Numerical Simulation 

In this Section, we discuss the numerical integration of the atomic and field equations 
derived given in the last Section. 

There are several possible ways to integrate the differential equation numerically. 
Three main categories of algorithm used are forward (explicit), backward (implicit), 
and mid-point (semi-implicit) methods [88]. The midpoint method is in a sense 
between the explicit and implicit methods, and we will use an algorithm of this type 
in the following. Let tm — tn + ^ for nth segment and iterate (m denotes mid point) 




(4.3) 
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until convergence is reached. Then step forward with = 2x{tm) — x{tn)- 

The forward difference method, which Euler or Runge-Kutta methods utihzes, is 
not guaranteed to converge in stochastic integrations [37]. There it is shown that 
the semi-imphcit method [95] is more robust in Stratonovich type SDE simulations 
[36] . More extensive studies of the stability and convergence of SDE can be found in 
the reference [96]. The Stratonovich type SDE equivalent to the Ito type equation 
(4.2), is 

j k 

+ ^Bij{t,x^t)dWf (Stratonovich), (4.4) 
j 

which has the same diffusion terms Bij, but with modified drift terms. This "cor- 
rection" term arises from the different definitions of stochastic integral in the Ito and 
Stratonovich calculus. 

At the end of Appendix B 3.3, we derive the Stratonovich SDE with the (un- 
derlined) "correction" terms noted above. We then have 19 classical variables in- 
cluding atomic populations, coherences, and two counter-propagating cascade fields. 
With 64 diffusion matrix elements and an associated 117 random numbers required 
to represent the instantaneous Langevin noises, we are ready to solve the equations 
numerically using the robust midpoint difference method. 

4.3.1 Shooting and secant method 

The problem we encounter here involves counter-propagating field equations in the 
space dimension and initial value type atomic equations in the time dimension. The 
initial value problem is addressed by the difference method discussed in the previous 
Section. 

The counter-propagating field equations have a boundary condition specified at 
each end of the medium. This is a two-point boundary value problem, and a numerical 
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initial guess 




z 



Figure 4.1: Schematic illustration of the principle of the shooting method for two- 
point boundary value problems. 
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approach to its solution, the shooting method [88], is illustrated in Figure 4.1. 

Consider the set of differential equations dXi(z)/dz = gi{z,X). A subset A of 
{Xi} satisfy boundary conditions at ^ = 0, and the complementary subset B satisfy 
boundary conditions at z — L. 

The shooting method augments the set A with a set of "guesses" A', so that AuA' 
enable the differential equations to be integrated as an initial value problem (from 
z — to z — L). The idea is that A' is the correct choice when the integrated 
values at z = L reproduce the true boundary conditions, set B, within a permissible 
tolerance. The set A' is updated to enable convergence of the output at ^ = L to 
the set B. 

The secant method that is used to update each element of A' takes two guesses 
Xi and X2 for each variable of A' and returns an updated value Xi, 

Xi^X2-f2^^^. (4.5) 
X2 - Xi 

where /i and /2 are the differences between the required values of that variable in 
set B and the numerically computed values assuming xi and X2 values at z = 0. 
This method is iterated until convergence to all values in B is obtained. The secant 
method is illustrated in Figure 4.2. 

4.3.2 Outline of the numerical solution 

We use Matlab to perform the numerical integrations. For simplicity, we label the 
atomic and field variables as and ej. The counter-propagating field {—z direction) 
variables are ei and 62 (signal fields) and 63 and 64 (idler fields) propagate in the +z 
direction. We set the local time r ^ t in the following description of the algorithm. 

We initialize 15 ai{z,t), 4 ei{z,t) in time t e (0,7") and space z e (0, L), and 
select 19 Gaussian random numbers ni{z, t). Set time and space grids with spacings 
At, Az respectively. For each realization among R statistical ensemble averages, we 
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Figure 4.2: Secant method. The root is bracketed by two initial guesses of Xi and 
X2 and an updated guess Xj is located at the intersection of two straight lines. 
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update the variables governed by the symbohc equations of motion, 



—ei = Pi{e,a,ne,), (4.6) 
d 

—ai^Ai{e,a,nai), (4.7) 

where and are in general the functions of variables that are denoted as vectors 
e and a. Each variable has its own stochastic source term as rig . or . . 

The algorithm proceeds by using the midpoint difference method for the evolutions 
in space and time and the shooting method for ej, 



eiizm,t) = ei{z,t) + —Pi[e{zm,t),a{z,t),nei{z,t)], 
ai(z,tm) = ai(z,t) + ^Ai[e(z,t),a(z,tm),nai(z,t)], 

where z^ — z + Az/2 and — t + At/2. The two guesses required in the secant 
method used in the shooting method are chosen as Xi — {ei(0, t), 62(0, t)} and X2 — 
{es{L,t),e4{L,t)}. 

Any normally-ordered quantity (Q) can be derived by ensemble averages that 
(Q) = J2f=i Qi/R where Qi is the result for each realization. Note that the update 
for field variables in space precedes the update for atomic variables, which takes into 
account that field variables evolve faster than atomic variables. The order should 
not matter when finer grids are used. 

4.3.3 Results for signal, idler intensities, and the second-order correlation 
function 

In this subsection, we present the second-order correlation function of signal-idler 
fields, and their intensity profiles. We define the intensities of signal and idler fields 

by 
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(4.8) 



respectively, and the second-order signal-idler correlation function 



G,,,{t,T) = {E;{t)Er{t + T)E^{t + T)Et{t)) 



(4.9) 



where r is the delay time of the idler field with respect a reference time t of the signal 
field. Since the correlation function is not stationary [64], we choose t as the time 
when Gg^i is at its maximum. 

We consider a cigar shaped ^^Rb ensemble of radius 0.25 mm and L = 3 mm. 
The operating conditions of the pump lasers are {fla, ^b, ^i, ^2) — (0.4, 1, 1, 0)703 
where is the peak value of a 50 ns square pulse, and fit, is the Rabi frequency of 
a continuous wave laser. The four atomic levels are chosen as (|0), |1), |2), |3)) = 
(|5Si/2,F==3), |5P3/2,F=4), |5P3/2,F=4), |4D5/2,F=5)). The natural decay rate for 
atomic transition |1) — >■ |0) or |3) — >■ |0) is 701 = 703 = 1/26 ns and they have a 
wavelength 780 nm. For atomic transition |2) — ?> |1) or |2) — > |3) is 712 = 732 = 
0.156703 [97] with a telecom wavelength 1.53/xm. The scale factor of the coupling 
constants for signal and idler transitions is gg/gi — 0.775. 

We have investigated four different atomic densities from a dilute ensemble with 
an optical density (opd) of 0.11 to a opd — 4.35. In Figure 4.3, 4.4, and 4.5, we 
take the atomic density p — 10^° cm"^ (opd = 2.18) for example, and the grid sizes 
for dimensionless time At = 4 and space Az = 0.0007 are chosen. The convergence 
of the grid spacings is fixed in practice by convergence to the signal intensity profile 
with an estimated relative error less than 0.5%. 

The temporal profiles of the exciting lasers are shown in the left panel of Figure 
4.3. The atomic density is chosen as p = 10^° cm~^, and the cooperation time Tc is 
0.35 ns. The right panel shows time evolution of atomic populations for levels |1), 
|2), and |3) a.t z = 0,L, that are spatially uniform. The populations are found by 
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Figure 4.3: Time-varying pump fields and time evolution of atomic populations. 
(Left) The first pump field fla (dotted-red) is a square pulse of duration 50 ns and 
Qb is continuous wave (dotted-blue). (Right) The time evolution of the real part 
of populations for three atomic levels cth = (ais) (dotted-red), (T22 = {0^12) (dotted- 
blue), (J33 = (ciii) (dotted-green) at z = 0,L, and almost vanishing imaginary parts 
for all three of them, indicate convergence of the ensemble averages. Note that these 
atomic populations are uniform ClS db function of z. 
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ensemble averaging the complex stochastic population variables. The imaginary parts 
of the ensemble averages tend to zero as the ensemble size is increased, and this is a 
useful indicator of convergence, see Appendix B.2 for a discussion. In this example, 
the ensemble size was 8x10^. The small rise after the pump pulse is turned off is 
due to the modulation caused by the pump pulse fib, which has a generalized Rabi 
frequency -\/A| + 40^. This influences also the intensity profiles and the correlation 
functions. 

In Figure 4.4, we show that counter-propagating signal {—z) and idler {+z) fields 
at the respective ends of the atomic ensemble. The plots show the real and imaginary 
parts of the observables, and both are normalized to the peak value of signal intensity. 
Note that the characteristic field strength in terms of natural decay rate of the idler 
transition (703) and dipole moment (tij) is {di/h)Ec w 86.8703. The fluctuation in the 
real idler field intensity at z — L and non-vanishing imaginary part indicates a slower 
convergence compared to the signal field that has an almost vanishing imaginary part. 
The slow convergence is a practical limitation of the method. 

In Figure 4.5 (a), we show a contour plot of the second-order correlation function 
Gs,iiis,ti) where ti > ts- In Figure 4.5 (b), a section is shown through ts ~ 75 ns 
where Gs,i is at its maximum. The approximately exponential decay of Gs,i is clearly 
superradiant consistent with the theory of Chapter 8 and the reference [16]. The 
non-vanishing imaginary part of Gs,i calculated by ensemble averaging is also shown 
in (b) and indicates a reasonable convergence after 8x10^ realizations. In Table 4.1, 
we display numerical parameters of our simulations for four different atomic densities. 
The number of dimensions in space and time is Mt x with grid sizes (At, A^) in 
terms of cooperation time (Tc), length (I/c). The superradiant time scale {Tf) is 
found by fitting Gg^i to an exponential function (e~*/^^), with 95% confidence range. 

In Figure 4.6, the characteristic time scale is plotted as a function of atomic den- 
sity and the factor Np,, and shows faster decay for optically denser atomic ensembles. 
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Figure 4.4: Temporal intensity profiles of counter-propagating signal and idler fields. 

(a) At 2; = 0, real (dotted-red) and imaginary (diamond-red) parts of signal intensity. 

(b) At z = L, real (dotted-blue) and imaginary (diamond-blue) parts of idler intensity. 
Both intensities are normalized by the peak value of signal intensity that is 7.56 x 10~^^ 
E^. Note that the idler fluctuations and its no n- vanishing imaginary part indicate a 
relatively slower convergence compared with the signal intensity. The ensemble size 
was 8x10^, and the atomic density p = 10^°cm~^. 
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Figure 4.5: Second-order correlation function Gs,i{ts,ti). The 2-D contour plot of 
the real part of Gs^i with a causal cut-off at tg = ti is shown in (a). The plot (b) 
gives a cross-section at tg = tm ^ 75 ns, which is normalized to the maximum of the 
real part (dotted-blue) of Gs^i- The imaginary part (diamond- red) of Gg^i is nearly 
vanishing, and the number of realizations is 8x10^ for p = 10^°cm~^. 
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Table 4.1: Numerical simulation parameters for different atomic densities p. Corre- 
sponding optical depth (opd), time and space grids (M^ x M^) with grid sizes (At, Az) 
in terms of cooperation time (Tc) and length (Lc), and the fitted characteristic time 
Tf for Gs,i (see text). 



p(cm ^) 


opd 


Mt x 


At(T,), 


Tc(ns), 
Lc(m) 


fitted Tj(ns) 
[95% confidence range] 


5xl0« 


0.11 


101 x 44 


0.9, 1.5x10-^ 


1.55, 0.46 


24.6 [24.2, 25.0] 


5x10^ 


1.09 


101 x 42 


2.8, 4.5x10"^ 


0.49, 0.15 


14.8 [14.4, 15.3] 


1x10^° 


2.18 


101 X 42 


4, 7x10-^ 


0.35, 0.10 


9.4 [9.2, 9.7] 


2x10^° 


4.35 


101 X 42 


5.5, 1x10-3 


0.24, 0.07 


5.0 [4.6, 5.5] 



We also plot the timescale Ti = /{N/jL + 1) (ns) that is derived from the theory 
of Chapter 3, in which p, is the geometrical constant for a cylindrical ensemble, Eq. 
(A. 12). The natural decay time 7^53^ = 26 ns corresponds to the D2 line of ^'^Rb. 
The error bar indicates the deviation due to the fitting range from the peak of Gs,i 
to approximately 25% and 5% of the peak value. The theory and simulations are in 
good qualitative agreement, approaching independent atom behavior at lower densi- 
ties. For larger opd atomic ensembles, larger statistical ensembles are necessary for 
numerical simulations to converge. The integration of 8x10^ realizations used in the 
case of p = 10^*^ cm^'^ consumes about 14 days with Matlab's parallel computing tool- 
box (function "parfor") with a Dell precision workstation T7400 (64-bit Quad-Core 
Intel Xeon processors). 

4-4 Conclusion 

We have derived c-number Langevin equations in the positive-P representation for the 
cascade signal-idler emission process in an atomic ensemble. The complete c-number 
Langevin noise correlations are derived and confirmed by an alternative theoretical 
method. The equations are solved numerically by a stable and convergent semi- 
implicit difference method, while the counter-propagating spatial evolution is solved 
by implementing the shooting method. 
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p (10^° cm"^) 



Figure 4.6: Characteristic timescales, Tf and Ti vs atomic density p and the super- 
radiant enhancement factor Nfi (/i = fi). Tf (dotted-blue) is the fitted characteristic 
timescale for Gs,i{ts = tm,ti = tm + t) where tm is chosen at its maximum, as in 
Figure 4.5. The errorbars indicate the fitting uncertainties. As a comparison, 
Ti='^Q^ /{Njji + 1) (dashed-black) is plotted where 7q3^ = 26 ns is the natural decay 
time of Dl line of *^Rb atom, and /i is the geometrical constant for a cylindrical 
atomic ensemble, as discussed in Chapter 3. The number of realizations is 4x10^ for 



P 



5 X 10^ 5 X 10^ cm-3 and 8x10^ for p = 10^°, 2 x 10^° cm" 
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We investigate four different atomic densities readily obtainable in a magneto- 
optical trap experiment. Signal and idler field intensities and their correlation func- 
tion are calculated by ensemble averages. Vanishing of the unphysical imaginary 
parts within some tolerance is used as a guide to convergence. We find an enhanced 
characteristic time scale for idler emission in the second-order correlation functions 
from a dense atomic ensemble, consistent with the superradiance timescales predicted 
by the analytical method in Chapter 3, and observed experimentally [16]. 
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CHAPTER V 



SPECTRAL ANALYSIS FOR 
CASCADE-EMISSION-BASED QUANTUM 
COMMUNICATION 

Cascade emission in alkali atoms is a source of telecommunication photons. In this 
Chapter, we investigate the DLCZ [6] scheme using the cascade emission from an 
atomic ensemble. 

5. 1 Introduction 

Long distance quantum communication based on atomic ensembles was proposed by 
Duan, Lukin, Cirac, and ZoUer [6]. This scheme involves Raman scattering of light by 
the atoms. The cascade transitions investigated in Chapter 3 and 4 provide a source 
of telecommunication wavelength photons. It is interesting to assess the cascade 
scheme in the DLCZ protocol given that it could potentially reduce transmission 
losses in a quantum telecommunication system. The DLCZ scheme is based on 
entanglement generation and swapping and quantum state transfer. 

In this Chapter, we first discuss entanglement generation and then investigate how 
frequency entanglement of the cascade photon pair influences entanglement swapping. 

5.2 DLCZ Scheme with Cascade Emission 

In the DLCZ protocol, a weak pump laser Raman scatters a single photon generating 
a quantum correlated spin excitation in the ensemble. By interfering the Raman 
photons generated from two separate atomic ensembles on a beam splitter (B.S.), the 
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DLCZ entangled state (|01) + |10))/-v^ [98] is prepared conditioned on one and only 
one click of the detectors after the B.S. Hence |0) and 1 1) represent the state of zero or 
one collective spin excitations stored in the hyperfine ground state coherences. This 
state originates from an indistinguishable photon paths. The error from multiple 
excitations can be made negligible if the pump laser is weak enough. 

As shown in Figure 5.1, we consider instead that one of the ensembles employ 
cascade emission. The idea is for cascade emission to generate a telecom photon (aj) 
for transmission in the optical fiber, and an infrared photon that interferes locally 
with the Raman photon generated in the A-type atomic ensemble. In this way 
interference of the infrared photons generate the entangled state, 

l*) = ^(|01)a,.+ |10)„,.), (5.1) 

similar to the conventional DLCZ entanglement generation scheme. Now, however, 
instead of a stored spin excitation, we generate a telecom photon. 

The entanglement swapping with the cascade emission may be implemented as 
shown in Figure 5.2, and will be discussed in detail in the next Section. The initial 
state is a tensor product of two state vectors generated locally at the sites A and B. 

I*) = (^A^|0) + VvU\i)t\i)f) (^A^|0) + Vv^\i)^\i)^)® 

{^/T^^\0) + Vm^|l)f |l)f ) ® {^/l^^\0) + V^|l)f |l)f ), (5.2) 

where (s, i) represent the signal and idler photons from the cascade emission, and (r, 
a) are Raman scattered photon and the collective spin excitation. Here rji and r]2 are 
efficiencies to generate cascade and Raman emission. Since rji and 772 <S 1, multiple 
atomic excitations or multi-photon generation can be excluded. 

5.3 Entanglement Swapping 

Consider the product state generated from A and B, Figure 5.2, 
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Figure 5.1: Entanglement generation in tlie DLCZ sclieme using the cascade and 
Raman transitions in two different atomic ensembles. Large white arrows represent 
laser pump excitations corresponding to the dashed lines in either cascade or Ra- 
man level structures. Here represents the emitted telecom photon. B.S. means 
beam splitter that is used to interfere the incoming photons measured by the photon 
detector D. The label A refers to the pair of ensembles for later reference. 
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Figure 5.2: Entanglement swapping of DLCZ sclieme using the cascade transition. 
The site A is described in detail in Figure 5.1 and equivalently for the site B. The 
telecom signal photons are sent from both sites and interfere by B.S. midway between 
with detectors represented by c\ and c^. Synchronous single clicks of the detectors 
from both sites (m| 2, n\ 2) and the midway detector (c| 2) generate the entangled state 
between lower atomic ensembles at sites A and B. The locally generated entanglement 
is swapped to distantly separated sites in this cascade-emission-based DLCZ protocol. 
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|10)., + |01). |10),. + |01). 
I*) - ( ^ )a ^ ( ^ )b 

= ^{\WW)asas + \WOl)asas + \OlW)asas + \QWl)asas), (5-3) 

where the subscript (a) represents a stored local atomic excitation, and (s) means a 
telecom photon propagating toward the B.S. in the middle. We can tell from this 
effective state that the first component (|1010)osas) contributes no telecom photons 
at all (two local excitations) and can be ruled out by measuring a "chck" at one of 
the middle detectors. The second and the third components have components of the 
entangled state of quantum swapping, and the fourth one is the source of error if the 
photodetector cannot resolve one from two photons. The error could be corrected 
by using a photon number resolving detector (PNRD) if other drawbacks like dark 
counts, photon losses during propagation, and detector inefficiency are not considered. 

Now we will formulate the entanglement swapping including the spectral effects 
discussed in Chapter 3. We ignore pump-phase offsets, assuming 50/50 B.S. and 
a symmetric set-up (771A = 771B = rji, 772A = V2B — V2) for simplicity. Expand the 
previous joint state, Eq. (5.2) and keep the terms up to the second order of 771,2 that 
can contribute to detection events (^12)^2); 



l*)e// 

= m(l - r72)|l),^|l)f |l)f |l)f + V2{1 - m)|l);^|l)^.|l)f 

VvMi-vi){i-V2)\i)t\i)f\i)^\i)Z + VvM^-vi){i-V2mf\i)i\i)f\i)f, 

(5.4) 

where the cascade emission state = f f {us, Ui)a\^{us)a\,{ui)\0)duJsdUi has the 

spectral distribution f{ujs,uji) as derived in Chapter 3. 

As shown in Figure 5.2, entanglement swapping protocol is fulfilled by measuring 
three clicks from the three pairs of the detectors respectively (?«| 2) 't-^; c^)- The 
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quantum efficiency of the detector is considered in the protocol, and we describe a 
model for quantum efficiency in Appendix C.l. We then use this model to describe 
photodetection events registered by non-resolving photon detectors (NRPD). Starting 
with the input density operator = |^)e//(^|, we derive the projected density 
operator, Eq. (C.16), conditioned on the three clicks of m\,h\, and c\ in Appendix 
C.2. We use the Schmidt decomposition of the projected density operator and assume 
a single mode for the Raman scattered photon. We find the un-normalized density 

(2) 

operator p^^j given in Eq. (C.16), 



.(2) _ 
Hout ~ 

-(2 - Vtmeff[^ + 2^ \- j |0)(0| + ^ VtVeff 



16 



(^SU0){0\SA + Sl,\0){0\SB)y (5.5) 

where rjt and rjeff are quantum efficiencies of the detectors at the telecom and infrared 
wavelengths respectively. The first term in Eq. (5.5) is the atomic vacuum state at 
sites A and B and contributes an error to the output density operator. The second 
term contains the components of the DLCZ entangled state. 

We can define the fidelity F, the success probability Ps of entanglement swapping 
of the entangled state |^)dlcz — {Sa + SI)\0)/V2, and the heralding probabifity Ph 
for the third chck as 



DLCZ 



^ — „ . ,(2) ^ ' W-'^J 



Ph = Pi + P2, Pi = P2 = (5.7) 
Ps^PixFi + P2X F2, Fi = F2 = F, (5.8) 



66 



where Pi^2 is the heralding probabihty of the single click from the midway detector 
(c| 2) as shown in Figure 5.2, and a trace (Tr) is taken over atomic degrees of freedom. 
The normalization factor J\f is calculated in Eq. (C.9) and is given by 



vK^-mf 2 , ^^2(1 -^2) 2 , ^i(i-^i) 



2 



We have used the following properties for the calculation of p^J^f and jV, 



dusduji\f{ujs,<^i)\'^ = 1, (5.10) 
where orthonormal relations in the mode functions are used, and 

dujsduj'^dujidujlf{uj'^,uj'i)f*{uj'^,uji)f{ujs,uji)f*{ujs,uj'i) ^^Aj. (5.11) 



J 

Note that the single mode spectral function for the Raman photon satisfies J du\^ 
1. 

The fidehty, heralding, and success probability become 



ry,(2-r/,)(l + E,A,')/2 + 2 



(5.12) 



77,77,(2 -r/i)(l + E, Ai)/2 + 277t 



Ph = (r-^^TAi (5-13) 



1 + A, /0,(u;,)0*(a;D**(c^.)**('^Dc?'^.c?'^^ 
^ (V^+IW ' ^'-''^ 



where ~ 1 and 77, = Vi/V2- 

The fidelity depends on a sum of square of Schmidt numbers in the denominator 
and the mode mismatch between the idler and Raman photons in the numerator. Let 
us assume that the Raman photon mode is engineered to be matched with the idler 
photon mode of the largest Schmidt number {(f)i{uji) in our case), which is required 
to have a larger fidelity (so is the success probability) compared to other modes. We 
may also compare the NRPD with the performance of PNRD in the midway detectors, 
then we have the fidelity, heralding, and success probability. 
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r?^r?t(2-r?0(l+E_,- A,")/2+2r)t 
(v^+l/^)2 
»jr>jt(l->?f)(l+EjA|)+2,7t 



r?r.(2-»?t)(l+E,A2)/2+2' 



(v^+l/v^)2 



7?.(l-»jt)(l+E,A|)+2' 



r;t(l+Ai) 



>?t(l+Ai) 



1+Ai 



1+Ai 



NRPD 



PRND 



, PRND 



, NRPD 



, PRND 



NRPD 



(5.17) 



(5.16) 



(5.15) 



When the relative efficiency is made arbitrarily small, the fidelity approaches 
(1 + Ai)/2 for both types of detectors. It reaches one if a pure cascade emission source 
is generated (von Neumann entropy E — and Ai = 1). When r}r — 1 with a pure 
source using NRPD with a perfect quantum efficiency, F — 2/3, Ph — ^/A, P5 = 1/2, 
which coincide with the results of the reference [99] (with perfect quantum efficiency). 

We discuss the frequency entanglement for various pulse widths and superradiant 
decay rates in Chapter 3.4. We find that for shorter driving pulses and smaller 
superradiant decay rates, the cascade emission source is less spectrally entangled. 
That means when 77^ is fixed, a shorter driving pulse heralds a higher fidelity DLCZ 
entangled state. 

In Figure 5.3, we numerically calculate the entropy and plot out the fidelity from 
Eq. (5.15), the heralding probability from Eq. (5.16), and the success probability 
from Eq. ( 5.17) as a function of the relative efficiency r]r- With a perfect detection 
efficiency [r] = 1), we find that at a smaller rjr, the less entangled source gives us a 
higher fidelity DLCZ entangled state but with a smaller success probability. Small 
generation probability for cascade emission (rjr < 1) reduces the error of NRPD from 
two telecom photons interference, but it reduces the successful entanglement swapping 
at the same time. 

The optimal success probability occurs by using the same excitation efficiency 
for both cascade and Raman configurations. For PNRD, the fidelity is higher than 
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Figure 5.3: Fidelity F, heralding Ph, and success Ps probabilities of entanglement 
swapping versus relative efficiency 77,. with perfect detection efficiency rjt = 1. Col- 
umn (a) NRPD and (b) PNRD. Solid-red, dashed-blue, and dotted-green curves 
correspond to the pulse width parameters r = (0.1,0.5,0.5) and superradiant factor 
NfL + 1 = (5, 5, 10) (see Chapter 3 and Appendix A). The von Neumann entropy is 
S = (0.684,2.041,2.886), respectively. 
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NRPD, and the heralding probabihty is the same independent of the degree of fre- 
quency space entanglement. The success probabilities for both types of detectors are 
equal. The advantage of PNRD shows up in the fidelity of quantum swapping. 

In Figure 5.4, we show that the measures improve monotonically with the quantum 
efficiency {rj = rjt) of the detector at telecom wavelength, with rjr = 0.5. The success 
probabilities for both types of detectors are the same and again the advantage of 
PNRD shows up in the fidelity. 

5.4 Polarization Maximally Entangled State (PME 
State) and Quantum Teleportation 

In Figure 5.5, we illustrate schematically a scheme for probabilistic PME state prepa- 
ration and quantum teleportation. Four ensembles (ABCD) are used to generate 
two entangled pairs of DLCZ entangled states, and another two ensembles (/i, I2) 
are used to prepare a quantum state to be teleported. 

With the conditional output density matrix from Eq. (C.16), wc proceed to con- 
struct the PME state \'^)pme = :^{SaSd+S^b^I)\0) where (C, D) represents another 
parallel entanglement connection setup. Figure 5.5 (a). This PME state is useful in 
entanglement-based communication schemes [6] , and we will here calculate its success 
probability. The normalized density matrix for the AB system is from Eq. (5.5) (let 
Vt = v), 

= ^I0)(0| + + S'a\0){0\Sa 

+ XiS],\0){0\Sa + XiSi\0){0\SB) , (5.18) 

where the largest Schmidt number (Ai) of mode overlap is chosen and a = r)r{2 — 
A parallel pair of entangled ensembles (C,D) is introduced, and the joint density 
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Figure 5.4: Fidelity F, heralding Ph, and success Ps probabilities of entanglement 
swapping versus telecom detector quantum efficiency 1] for the case of (a) NRPD and 
(b) PNRD. Solid-red, dashed-blue, and dotted-green curves correspond to the same 
parameters used in Figure 5.3. 
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Figure 5.5: PME projection (a) and quantum teleportation (b) in the DLCZ scheme. 
Four atomic ensembles (A,B,C,D) are used to generate two DLCZ entangled states 
at (A,B) and (C,D). PME state is projected probabilistically conditioned on four 
possible detection events of (D^ or dIj) and or D\j) in (a). In the quantum 
teleportation protocol (b), another two ensembles (Ii,l2) are used to prepare a quan- 
tum state that is teleported to atomic ensembles B and D conditioned on four possible 
detection events of {Dj^ or Da) and {Di^ or Dc)- 
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operator is Poui,n ® Pout,n ■ The latter expression is developed mathematically in 
Appendix C.3. 

With projection of the PME state, we have the post measurement success proba- 
bility [a chck from each side; the side of (A or C) and (B or D)], 



42),AB „ -(2).CD|,-r,\ 
S,PME - {^\Pout,n ^ Pout,n \^)PME, 

4(1 + A?) 



(5.19) 



For rjr <^ 1, Ps,PME reaches the maximum of 1/2 when a pure source (Ai = 1) is 
used. 

For an arbitrary quantum state transfer to long distance, quantum teleportation 
scheme may be used. Another two ensembles (Iijh) are introduced [6], and the 
quantum state can be described by |^) = {d^Sl^ + diS\^)\0) with \do\'^ + \di\'^ ~ 1. 
The joint density matrix for quantum teleportation is 

PQT = {doSl + d,Si)\0){0\{d^,Sj, + dlSj,) ® pS,f ® pSJ^. (5.20) 

Atomic ensembles (A,B) in parallel with (C,D) provide a scheme for PME state 
preparation. Retrieve the quantum state [ensemble (Ji, I2)] into photons and interfere 
them at B.S., respectively, with photons from A and C. We have the teleported 
quantum state at B and D conditioned on the single chck of {Dj^ or Da) and {Dj^ or 
Dc). 

Consider single detection events at Dj^ and Dj^ as an example. With the NRPD 

measurement operators M/i, 72 = (4i-|0)i3i(0|)«)|0)D^(0|®(/i,2-|0)D2(0|)(g)|0)Dc(0| 
(we use Di, D2 for Dj^, Dj^), the density matrix after the measurement becomes 
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Pi = Tr{pQT,effMi,,i^) = 

a + 2 4 /Mop At i^v /^i ^ Ml 



\0)abcd{0\ + 7^^(^'^l|0)(0|'^B + ^^L|0)(0|^,,+ 



2(a + 4)2'"'"""^^"' (a + 4) 
^4|0)(0|^^ + ^^L|0)(0|^^), (5.21) 

where PQT,eff is calculated in Eq. (C.18), and the trace is taken over the electromag- 
netic field degrees of freedom. 

For a successful transfer of the quantum state |$) = (do'S'^-|-(ii5'|))|0), the fidelity 
Fi = ($|pi|<l>)/Tr(pi), and the heralding probability is Pi = Tr(/}i), with the trace 
over all atomic degrees of freedom. Except for the detection event we consider here, 
there are three other detection events including {Da, Dc), (-D/j, Dc) and {Da, Dj^). 
The teleported state from the detection events (-D/^, Dq) and {Dj^, Da) requires a 
TT rotation correction on the relative phase (do do, di —di). 

The fidelity and heralding probabilities conditioned on the other three pairs of 
clicks are the same as Fi and Pi respectively, so the success probability is 

4 
i 

-J^^^2i^ + i'^>^l-'^)\do\'\din (5.22) 

where F is the fidelity of entanglement swapping for NRPD, Eq. (5.15). For PNRD, 
the success probability for quantum teleportation is unchanged. 

The success probability for quantum teleportation depends on the probability 
amplitude of the quantum state and the fidelity F of the entanglement swapping. In 
Figure 5.6, for r]r — 0.5 and rjt — 1, we can see in the region Mo| ~ 0.3 ~ 0.9, higher 
success probability requires a less entangled cascade emission source. Outside this 
region, it prefers a more entangled source. When a pure source is used (Ai = 1) and 
let rjr <^ 1, rjt — 1, we can achieve the maximum of the success probability Ps,qt — \ 
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Figure 5.6: Success probability of quantum teleportation as a function of the prob- 
ability amplitude of teleported quantum state with rjr = 0.5 and a perfect detector 
efficiency rjt = 1. Solid-red, dashed-blue, and dotted-green curves correspond to the 
same parameters used in Figure 5.3. 

when F = 1, which is also achieved in the traditional DLCZ scheme with perfect 
quantum efficiencies [99] . 

5.5 Conclusion 

We have described probabilistic protocols for the DLCZ scheme implementing the 
cascade emission source. We characterize the spectral properties of the cascade emis- 
sion by Schmidt mode analysis and investigate the fidelity and success probability 
of the protocols using photon resolving and non-resolving photon detectors. The 
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success probability is independent of the detector type, but photon number resolving 
detection improves the fidelity. 

The performance of the protocol also depends on the ratio of efficiencies in gen- 
erating the cascade and Raman photons. The success probability is optimized for 
equal efficiencies while the fidelity is higher when the ratio is smaller than one for 
non-resolving photon detectors. 

The frequency space entanglement of telecom photons produced in cascade emis- 
sion deteriorates the performance of DLCZ protocols. The harmful effect can be 
diminished by using shorter pump pulses to generate the cascade emission. A state 
dependent success probability of quantum teleportation was calculated, and in some 
cases a more highly frequency entangled cascade emission source teleports more suc- 
cessfully. An improved performance could be achieved if the error source (vacuum 
part) were removed. This could be done by entanglement purification [3] at the stage 
of entanglement swapping and then using the purified source to teleport the quantum 
state. 
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CHAPTER VI 



EFFICIENCY OF LIGHT-FREQUENCY 
CONVERSION IN AN ATOMIC ENSEMBLE 

In this Chapter ^, the efficiency of frequency up and down conversion of light in an 
atomic ensemble, with a diamond level configuration, is analyzed theoretically. The 
conditions of pump field intensities and detunings required to maximize the conversion 
as a function of optical thickness of the ensemble are determined. The influence of 
the probe pulse duration on the conversion efficiency is investigated by the numeric 
solution of the Maxweli-Bloch equations. The set of equations are similar to those 
in Chapter 4, but a c-number version of the interaction is considered here. The 
properties of absorption and dispersion of fields are extracted from the steady state 
solutions to demonstrate the parametric coupling between the fields. We will show 
that, in calculating conversion efficiency, a quantum version of the equation including 
Langevin noises is equivalent to the c-number one. Frequency conversion provides 
the bridge for transmitted qubit (telecommunication wavelength) and local quantum 
memory (near-infrared light) , in which a large scale quantum communication can be 
fulfilled. 

In Section II, we discuss the four-wave mixing process and present solutions for 
the up- and down- converted fields. The dressed state picture is used as a guide to 
understand the characteristic features of the absorption and signal-idler field coupling. 
In Section III, we present the results of an optimization in conversion efficiency as a 
function of the optical depth of the atomic ensemble. In Section IV, we investigate 

^This Chapter is based on reference [100]. 
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the effects of a finite pulse duration by numerically integrating the Maxwell-Bloch 
equations. Section V demonstrates the results of Langevin noise correlations and 
we conclude in Section VI. The derivations of the Maxwell-Bloch and parametric 
equations are relegated to the Appendix D. 

6. 1 Introduction 

The frequency conversion of light fields has been an important theme in optical physics 
for around half a century. In quantum information physics the conversion of single 
photons to and from the telecom wavelength band is a topic of more recent vintage, 
and is motivated by the desire to minimize optical fiber transmission losses when 
distributing entangled states over distant quantum memory elements in a quantum 
repeater [4]. 

An associated technical problem is that telecom light is not readily stored in 
ground level atomic memory coherences. Retrieval processes in atomic ensembles, 
for example using electromagnetically induced transparency [101], or more specifi- 
cally the dark-polariton mechanism [102, 103], generate shorter wavelength radiation 
correlated to the stored atomic excitation by Raman scattering. Such radiation, op- 
tically resonant to the ground level of typical atoms and ions, has been retrieved in 
numerous experiments [7, 8, 10, 9, 11, 13, 14, 15, 104, 105]. An important advance 
would involve generation of atomic memory coherences quantum-correlated with tele- 
com wavelength radiation, thereby minimizing transmission losses over long distances. 
Recently there has been a breakthrough in this direction using a pair of cold, non- 
degenerate rubidium gas samples [25]. The stored excitation is correlated with an 
infra-red field (idler) in one gas sample, and the idler is then frequency converted to 
a telecom wavelength signal field in the other ensemble. The frequency conversion 
mechanism involves the diamond configuration of atomic levels shown in Figure 6.1. 

In a probabilistic protocol it is important to maximize all efficiencies, e.g., fiber 
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Figure 6.1: The diamond configuration of atomic system for conversion scheme. 
Two pump lasers (double line) with Rabi frequencies fla, and propagated probe 
fields (single line) E^, interact with the atomic medium. Various detunings are 
defined in the Appendix D, and the atomic levels used in the experiment [25] are 
(|0), |1), |2), |3)) = (|5Si/2, F = 1), I5P3/2, F = 2), |6Si/2, F = 1), |5Pi/2, F = 2)). 

transmission, single-photon detection, and quantum memory lifetime [106]. In the 
present work we investigate the efficiency of frequency up- and down- conversion in 
the diamond atomic configuration [16, 107], as a function of the ensemble's optical 
thickness, and the intensity and detuning of the pump fields involved in the near- 
resonant, four- wave mixing process. 

6.2 Theory 

We consider a cold and cigar-shaped ^^Rb atomic ensemble with co-propagating light 
fields similar to the experimental setup in the reference [25]. 

The conversion scheme shown in Figure 6.1 involves two pump lasers with fre- 
quencies Ua and Uh, respectively; their Rabi frequencies are given by Qa and Qb. Two 
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weak probe fields, signal and idler, with frequency Ug and Ui, respectively, propa- 
gate through the optically thick atomic medium. Unlike the cascade driving scheme, 
where two-photon excitation generates a photon pair spontaneously [16], pump laser 
b experiences a transparent medium if both the signal and idler fields are in the vac- 
uum state. With an incident signal field, four-wave mixing with the pumps generates 
an up-converted idler field, while an incident idler field generates a down-converted 
signal. 

The Maxwell-Bloch equations for the interacting system of light and four light 
fields is derived in the Appendix D. By linearizing the equations with respect to 
the signal and idler field amplitudes, and adiabatically eliminating the atoms, one 
arrives at coupled parametric equations for the signal and idler fields. We discuss 
their solution in this section, and leave numerical solutions of the Maxwell Bloch 
equations to Section IV. 

The calculation of conversion efficiencies can also be carried out with the quan- 
tized Heisenberg-Langevin version of the coupled parametric equations, which we will 
show in Section V. The resulting conversion efficiencies are identical to the semiclas- 
sical treatment; the additional quantum noise contributions vanish as the |2) — )■ |3) 
transition driven by pump laser b has vanishing populations and atomic coherence. A 
similar simplification occurs in the calculation of the storage efficiency of spin waves 
in a system of atoms in the A configuration [108, 109]. 

The co-moving propagation equation for c-number signal and idler fields (respec- 
tively, E'^ and E^) under energy conservation {Aou — oua + oug — oub — oui — 0) and 
four- wave mixing conditions {Ak — ka — kg + kh — ki = 0) are 



d_ 

dz 

d_ 

dz 



Et = ^sEt + KgEt 




KiEj + aiEj. 



(6.1) 



80 



The coupled equations are similar to those found for the double A system [110, 
111]. The self-coupling coefficients (3s, oti and parametric coefficients defined 
in Appendix D. The set of equations can be simplified as 



d_ 
dz 



x{z) = Ax 



(6.2) 
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The equations are solved by considering a similarity transformation S that A = 
S~^AS is diagonahzed and y — S~^x such that 

d 



(6.4) 
(6.5) 



where y{zo) is the boundary condition. With the known boundary condition Xi{0) 
and ^2(0) where we choose the input face of propagation as Zq — 0, we have 



x{z) = Se^'S-^x{0). 



(6.6) 



And the diagonahzed and transformation matrix are 
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where w = \/ q^ + KsKi, and q = {—ai + /3s) f^,. 
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The solution of fields from down conversion is 
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Similarly, the solution of fields from up conversion is 
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We define the down conversion efficiency t]^ and transmission of input idler field 
Td as 
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(6.12) 
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Similarly the up conversion efficiency r]^ and transmission of input signal field Tu 
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(6.14) 
(6.15) 



2w{w + q) 

The above is the central result of this section. The up and down conversion 
efficiencies differ only in the parametric coupling coefficients and Kg- In the 

strong parametric coupling regime where IkjUk^I ^ |Q;j|, 1/3^1, the coefficients can 
be simplified to — \ sinh(y'ACsACjL), 77^ ~ sinh (^Kg^jL) and Tu = ~ 
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Figure 6.2: Dressed-state picture from tlie perspective of tlie probe idler transition 
between atomic levels |0) and |3). Two strong fields Qa, shift the levels with energy 
AEa^b and wavy lines represent the idler field resonances. 

cosh.{y/ KgHiL). Under the further assumptions = /3s = and Kj, Ks are pure imag- 
inary, we find r|^^ = r]^ = siri^llm^KsL)] and T^^ = = cos'^[lm{KsL)]. This result 
was recently derived by Gogyan using a dressed state approach [112], in the case of 
resonant pump fields Ai = A;, = [113]. In this ideal limit there is a conservation 
condition rj^+T^^ = r^d+^d = 1- The parametric coupling coefficients are not identical, 
but in the regime of strong coupling they approach each other. As noted by Gogyan, 
when the pump-a intensity is large (Qa » |Ai|,7o3) the |0) ^ |1) is saturated, the 
atomic coherence is negligible and Kg ^ Hi oc 5"oo,s(^^^ + ^^j^)- Alternatively, in the 
limit fib » ^a,732 and |Ai| >> 701 the atomic coherence of |0) ^ |1) dominates 
and once again ~ oc —^^Ifr^^- Note that this scheme is also similar to the 
frequency conversion in nonlinear materials [114]. 
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The ac-Stark splitting induced by the pump lasers shifts the resonant absorp- 
tion condition for the idler and signal fields. The idler and signal experience reso- 
nant absorption at the transition frequency of the dressed atom. The corresponding 
transitions for the idler are shown in Figure 6.2. The bare states are shifted by 
AEa = Ai ± ^Af + 4QI /2 and AEb = A^ ± ^A^ + AQl /2, respectively Note 
that our Rabi frequencies are smaller by a factor 2 than the standard definitions to 
avoid a plethora of prefactors in the equations of the Appendix. 

For resonant pump fields, AEa,b = ±f^a,6 • The idler transition resonances are 
at Auji — —{^a + ^b), — — ^b\ , — ^b\ , (^a + ^b) and thesc delineate three 
windows separated by these four absorption peaks. For > ilj, the centers of these 
windows are at — fla, 0, and fla, respectively. Choosing the idler detuning AoUi = iQ^ 
as in Ref. [113], the idler interacts with the atomic medium at the center of the left 
or right window. 

As an example of the strong coupling windows created by intense pump lasers, 
we show in Figure 6.3 the self and cross coupling coefficients for the signal and idler 
fields as a function of the idler frequency. Note that the corresponding frequency of 
signal field is determined by Aus = Aui — Ai + A;,. The dimensionless quantities 
aiL, and KsL are shown under the conditions of maximum conversion efficiency 
to be discussed in the next section. We choose the optical depth (opd) paL — 150 
where p is the number density, a = 3A^/ (Att) the resonant absorption cross-section, 
and L the atomic ensemble length in the propagation direction. Three parametric 
coupling windows are separated by two strong absorption peaks on the left and two 
relatively weak ones on the right. The imaginary part of the self-coupling coefficients 
are seen to vanish in each window at a certain point, while the real parts are small 
away from resonances. At the same time the cross-coupling coefficients have a large 
imaginary part. The positive gradient of ^(/^^L) and '^{aiL) inside the windows is 
indicative of normal dispersion. 
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Figure 6.3: Self-coupling coefficients (3^, ai and cross-coupling coefficient Hg- Dimen- 
sionless quantities (a) (3sL, (b) aiL and (c) k^L with real (solid blue) and imaginary 
(dashed red) parts are plotted as a dependence of idler detuning Aw, [same label 
in (b)] showing four absorption peaks to construct three parametric coupling win- 
dows. A black dashed-dot line of the constant 7r/2 is added in (c) to demonstrate 
the crossover with 9(ksL) indicating the ideal conversion efficiency condition in the 
left window. The parameters we use are (f2a, f^b, Ai, A;,, Acjj) = (33, 20, 39, 2, 
— 21)703 for optical depth paL = 150 with L = 6mm. Various natural decay rates are 
7o3 = l/27.7ns, 701 = 1/26. 24ns, 712 = 703/2.76, and 732 = 7o3/5.38 [97]. 
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Figure 6.4: Down conversion efficiency r^d vs optical depth (opd) from 1 to 600. 
Each dotted point is the maximum for five variational parameters Qa, ^b, ^i, ^b, 
and Aui. 

6.3 Optimal Conversion Efficiency 

It is important to ascertain the parameters that allow maximum efficiency of conver- 
sion due its potential in practical quantum information processing. In principle we 
need to search the three parametric coupling windows to find the optimum conditions 
for an atomic ensemble of a given optical thickness. 

In the previous section we have discussed how three parametric coupling windows 
appear for some particular values of pump laser parameters. In the search for the 
maximal conversion efficiency, five parameters fi^, ^b, and Aui are varied to 

maximize the conversion efficiency for a fixed optical depth of atomic ensemble, using 
functional optimization. 
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The optical depth puL appears through the dependence on atomic number N in 
the Arecchi-Courtens cooperation time [115] 

In Figure 6.4, we show the maximum of down conversion efficiency using Eq. (6.12) 
for different optical depths from 1 to 600. The maximum is found by varying five pa- 
rameters mentioned above and the conversion efficiency reaches 100% asymptotically 
when the optical depth becomes larger. In the strong parametric coupling regime 
as we discussed in the previous section, r]^ ~ sin^[Im(«;5L)] and it has a maximum 
when \m{KsL) = |, see Figure 6.3. Since Im(KsL) is proportional to optical depth 
and inversely proportional to the Rabi frequencies of the driving lasers, an order of 
magnitude estimate of the optical depth necessary for near unit conversion efficiency 

is 0pd~ fl^a,6/703 >> 1- 

The behavior of the cross-coupling coefficient Im(KsL) as a function of idler de- 
tuning indicates where large conversion is to be found, as a comparison with Figure 
6.5 shows. The maximum efficiency of about 0.92 is located in the left parametric 
coupling window at the intersection of lm{KsL) and |. Inside the windows the trade- 
off between conversion and transmission is clear. In the region where absorption is 
large, on the sides of the window (especially for the left window), the efficiency and 
the transmission are both low although the valley in conversion efficiency corresponds 
to a peak in transmission as expected in parametric coupling. The transmission ap- 
proaches unity when the incident idler field is far off-resonance. 

We note that the symmetry (Ai, A;,, Aa;i) —{Ai,Af„AuJi) gives degenerate 
optimal conversion conditions. 

Moreover, for the region where absorption is large on the sides of the window 
(especially for the left window), the efficiency and the transmission are both low 
but the valley of efficiency corresponds to a peak for transmission indicating the 
feature of parametric coupling. The plateau for efficiency in absorption region can 
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Figure 6.5: Conversion efficiency rj^, rju and transmission Td vs Acuj for opd=150. r/d 
and r/u are indistinguishable and shown in solid red line, and Td is in dashed blue line. 
High transmission efficiency corresponds to low conversion efficiency indicating the 
approximate conservation condition within each parametric coupling window. The 
maximum conversion efficiency is found in the left window at around Aui = — 2O703 
and other relevant parameters are the same as in Figure 6.3. 
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be estimated as r^d ~ |^| where ctj^^ fa Ksi^i- Based on the study of finding the 
maximum efficiency for frequency conversion, wc will investigate the situation when 
the input is a pulse and numerical integration of full equation of motion is required 
in the next section. 

6.4 Pulse Conversion: Solution of the Maxwell- 
Bloch Equations 

The effect of finite-duration input probe pulses, which are often employed in practice, 
can be assessed by numerically solving the Maxwell-Bloch equations for the coupled 
atoms-fields system. The characteristic scales of time and length are given by the 
Arecchi-Courtens time and — cT^, respectively, which are inversely proportional 
to the square root of the opd. The cooperative electric field is the product of the 
atomic number and the idler electric field per photon, i.e., Ec — phhJij (2eo). 

Scaling the space, time, electric field amplitude, various detunings, and natural 
decay rates accordingly, indicated by tildes, the Maxwell-Bloch equations of Eqs. 
(D.4,D.6,D.7) under energy conservation (Ao; — uJa^ojg—^h — 'jJi — 0) and four- wave 
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mixing conditions {Ak — ka — kg + kf, — ki — 0) become 

9 ^ r^A 7ol^~ 



= (iAi - -^)aoi + i^ai^oo - ^u) + i^02Eg - ia\^E^, 



V . 2 

— Cri2 = (iAus - ^ ^^ )g"l2 - i^*ad-02 + ii^u - a22)E^ + i^bd-13, 
^5-02 = (iA2 - y)5-02 - iai2^a + i^oiEg + iaos^b - i^32El, 

^5-22 = -72 ^22 + io^i2Et - i^i2E~ + i0-hO^l2 - i^b^32, 

^C^33 = -703C^33 + 732C^22 " ^f^foC^L + ^^65-32 + l^lzE^ - iaQ^E^ ^ 

^cri3 = {i/\uji - lAi - ^^°'^ )ai3 - «n*ao3 - ^0-^2^5^ + ^^6^12 + ial^El, 

— ^03 = (^AWj - — )5-03 - if^a^l3 + «^^6C^02 + i(^00 " C^33)-Ej^, 

^c^32 = (-^^6 - - icri3E7 + ir^ft(a22 - c^33) + i^^E^ , (6.16) 

and 

where z = z/Lc, f = r/Tc, Qa,6 = ^a,bTc, E^^i = E^jEc, and [g'slVlfl'iP is a factor of 
unit transformation from signal to idler field strength. Natural life time [97] for signal 
and idler transitions is used to calculate the ratio of coupling strength gs/gi — 1.035. 
The above equations were integrated with a semi-implicit finite difference method 
[95]. The midpoint integration method is stable and has high accuracy without 
sacrificing memory for finer grids [88] . The algorithm has been tested by comparing 
with the parametric equations' solutions in appropriate limits, and these solutions 
arc recovered when fine enough grids are employed. 

To illustrate the infiuence of finite pump pulse duration, we compute the down 
conversion efficiency 

_ jmz = L,r)Pdr 
* - /|E.-(. = 0,r)|Mr- (^-'^^ 

In Figure 6.6, we show the computed values of r)d for two different input idler 
pulse durations. We fix the opd=150 and use the near optimum parameters (^2^, 
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Figure 6.6: Time-varying pump fields of Rabi frequencies fla,b{t) and down- 
converted signal intensity {\E^{t, z = L)\'^) from an input idler pulse {\E^{t, z = 0)p). 
Here we let t = r that is the delayed time in co-moving frame. Pump-b (dotted green) 
is a continuous wave and pump-a (dashed black) is a square pulse long enough to en- 
close input idler pulse with (a) 100 ns and (b) 15 ns (dashed-dot blue). Output 
signal intensity (sohd red) at the end of atomic ensemble z = L is oscillatory due to 
the pump fields. The square pulse in rising region (tr — ^ < t < tr + ^) has the form 
of |[1 + sin(^^^^)] that in (a) (t^, t,) = (10,10)ns for pump-a and (t^, t,) = (20,20)ns 
for input idler; (b) (t^, is)=(10,5)ns for pump-a and (t^, is) = (15,10)ns for input idler 
where tr is the rising time indicating the center of rising period tg- Note that the 
falling region of square pulse is symmetric to the rising one. 
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Qb, Ai, Afe, Auji) — (33, 20, 39, 2, —21)703 determined from the coupled parametric 
equations. The temporal shape of the pump laser intensities is also shown. Pump 
laser b is taken to be continuous wave, while pump a is a square pulse with duration 
large enough to completely overlap the input idler pulse. To compare with the 
steady state solutions, we choose the Rabi frequency of idler as O.I703, which is small 
compared to those of the pumps. We find that the conversion efficiency is reduced for 
shorter idler pulse inputs. A 100 ns idler pulse is long enough that it has a almost the 
same maximum conversion efficiency of 0.92 as in Figure 6.4 for opd=150. While for 
the shorter idler pulse of 15 ns, the signal develops significant temporal modulation, 
and this reduces the conversion efficiency, although it is still quite appreciable. The 
modulation frequency is at the generalized Rabi frequency of pump-a + AQ^. 

We note the characteristic time and space scales of the calculations are T(. = 0.086 ns 
and Lc — 26 mm for a moderate atomic density p — 1.7 x lO^^cm"^ and L — 6 mm. 
The grid size for dimensionless time At = 0.5 and space Az — 0.001 were chosen 
for both 100 and 15 ns idler pulse durations, and the convergence is reached with an 
estimated relative error less than 1%. 

Moreover, we show in Figure 6.7 of three dimensional plots of signal \E^{z, t) p and 
idler intensities \E^{z, t)\'^. A 100 ns input idler pulse is demonstrated in time-space 
propagation which is converted to signal pulse at the output surface of the ensemble. 

6.5 Discussion of Quantum Fluctuations 

In this Section, we derive quantized Heisenberg-Langevin equations by adding corre- 
sponding Langevin noises to the coupled equations. Similar to the results of Appendix 
D, we have 



d_ 

dz 
d_ 

dz 



Et = ^Et + KsEi + /, 



Ef ^ KiEt + aiEf + fi 



(6.20) 



(6.19) 
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Figure 6.7: Three-dimensional line plots of converted signal and input idler inten- 
sities in t (ns) and L (mm). Here we let t = r that is the delayed time in co-moving 
frame, and the parameters are the same as in Figure 6.6 (a). 
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where signal and idler fields (respectively, E'^ and E^) are now quantized, and 
Langevin noises and /j) in linearized equations are 



f - '^^^ \(T 
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J 13 
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(6.22) 



where various atomic and field Langevin noises J^y are associated with coupled equa- 
tions of atomic operators and field operators E^^ when y — s,i. 

Note that various normal correlation functions of quantum Langevin noises have 
(^J^j{t, z)Tj{t', z')^ — ^S{t — t')5{z — z')Dij in continuous hmit. If ensemble average 
is taken over the above field equations, with the property of Langevin noises that 
(j-l{t^ z)^ = (Pi{t' ^ z')^ = 0, the field equations are reduced to c- number ones where 
fiuctuations due to Langevin noises do not matter. For calculation of normally- 
ordered operators, semi-classical approximation is valid even in quantum regime. We 
demonstrate in the following and include Langevin noises in derivation of solutions 
of field operators. 

The set of equations can be written as 

d 



dz 



x{z) ^ Ax + f 



(6.23) 



where x{z) 
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fs{z) 



Consider a similarity transformation 5" that 5" ^^5" = A and y — S ^x, then we 
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have 
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y^Ay + S-'f, 



y{z) = e^^'-^'^yizo) + / dz'e^^'-''^S-^f{z') 

J zo 

With the boundary condition Xi{0) and 0:2(0), we have 



x{z) = Se^'S-^x{0) + / dz'Se^^'-''^S-^f{z'). 
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So we have correspondingly 
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where w = q^ + KgKi, and q = (— ctj + /3s)/2. 

The solutions of fields from down conversion for example are 
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where rja and are conversion efficiency and transmission for down conversion as 
derived in Section II, arg represents tfie argument for complex numbers, and we see 
tliat extra terms involve Langevin noises, i and i are defined in the below. 



— ^ .[{w + g)2e[K+/3»)/2+-](^--') + t,^^,e^(-'+M/^-^](^--\ (6.29) 

2w{w + q) 



^. ^ ^U{ai+M/2+w]iL-z') _ ^[(a,+M/2-w]{L-z')l /g_3Q^ 

= _^!i[e[K+/3.)/2+«'](i-^') _ g[(ai+/3.)/2-«;](L-z')]^ (6 31) 

2w[w + q) 

The down conversion efficiency 77down is defined as 
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_ {E:{L)Et{L)) 

"^'^ {EriO)Em) 
_ 1 

(^f^ dm{z')n{z')+U^')fl{z')] dz'%{z'')h{z'')+U^'')h{z'')]^ 

(6.33) 

in which expectation value of normally-ordered operators involves contributions from 
semi-classical treatment and normally-ordered noise correlation functions. The rele- 
vant normally-ordered quantum diffusion coefficients Dij from Einstein's relation are 
(note that Dij = Z)Jj 



(i) -Di2,i2 = 7oi {^22) ~ 7oi^22,s = 0; 
-Di2,o3 = -^12,02 = 0; 

-Di2,i3 = 701 (c^23) ~ 7oi^23,s = 0; 

Di2,s = Di2,i = 0; (6.34) 

(ii) £)o3,03 = 732 (C^22) ~ 732^22,s = 0; 

-Do3,02 = -Do3,i3 = 0; 

^03,. = Do3,i = 0; (6.35) 
(iii) I)o2,o2 = -Do2,i3 = 0; 

^02,. = Do2,i = 0; (6.36) 
(iV) Z)i3,i3 = 7oi (0-33) + 732 (5-22) ~ 7oi5-33,^ + 7325-22,5 = 0; 

^13,. = Di3,i = 0; (6.37) 

(V) Ds,s = Ds, = 0; (6.38) 

(Vi) = 0; (6.39) 
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where we have approximated various nonvanishing quantum diffusion coefficients by 
zeroth order properties of atomic operators (the steady state solutions). The above 
normally-ordered correlation functions give zero contributions in the linearized equa- 
tions of motion, so c-number Langevin equation is sufficient to derive the conversion 
efficiency. The normally-ordered noise correlations are zero because the population 
(^22,s' ^33,s) coherence {&23,s) properties are zero for the atomic level driven by 
pump-b. The linearized field equations in diamond structure have similar noise prop- 
erties to A system in which most atoms are on the ground state, and Langevin noise 
can be neglected if normally-ordered quantities, say storage efficiency, are considered 
[108, 109]. 

6. 6 Conclusion 

We have studied light frequency conversion in an atomic ensemble with a diamond 
configuration of atomic levels such as ^^Rb. The motivation stems from the need to 
efficiently convert light resonant with ground state transitions (storable in the sense of 
quantum memories) to and from the telecom wavelength band for low-loss quantum 
network communication. The optically thick atomic sample is driven by two strong 
co-propagating pump fields, and a probe idler or signal field depending on whether 
we consider down- or up-conversion. Parametric equations for the probe fields are 
derived and used to compute conversion efficiencies. They can be understood by 
dressed-state picture where we can visualize four absorption lines due to two strong 
pump lasers and thus three parametric coupling windows are created. There are two 
major contributions to the conversion efficiency, which are related to atomic popula- 
tions and coherences in the lower arm of the diamond level driven by laser pump-a. 
When this transition is saturated by a large pump Rabi frequency or when the co- 
herence dominates due to a large pump-b Rabi frequency in the upper transition, the 
cross-coupling coefficients and hence the conversion efficiencies are equal. 
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By performing a global parameter search we find conditions of pump Rabi fre- 
quencies, detunings, and signal/idler input frequency to maximize the conversion 
efficiency as a function of optical depth of the ensemble. Only in the limit of very 
large optical depth does the maximum efficiency approach the ideal strong coupling 
result [113]. Under conditions routinely obtained in cold, non-degenerate rubidium 
gas, with opd ~ 100 — 200, optimal conversion efficiencies of the order 80% to 90% are 
predicted. Numerical solution of the Maxwell-Bloch equations confirms the solution 
of the parametric equations in the limit of long pulse duration, and indicates that for 
shorter pulses, pump pulse induced modulation may reduce the conversion efficiency. 
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CHAPTER VII 
CONCLUSION 

We provide a theoretical study of light-matter interactions in cascade and diamond 
type atomic ensembles. A correlated two-photon (telecom signal-infrared idler) state 
vector is derived in the long time hmit within the adiabatic approximation. The 
second-order correlation function is calculated, and shows a superradiant time scale 
in the infrared idler emission. The entanglement in frequency space for such a two- 
photon state is analyzed by Schmidt decomposition. We are able to derive the mode 
functions and investigate the influence of pump pulse duration and superradiant decay 
rate that depends on optical density and ensemble geometry. 

To investigate multiple atomic excitations on the correlated emission from the 
atomic cascade transitions, we use the coherent state positive-P representation and 
derive an equivalent Ito type stochastic differential equation (SDE). The equations 
are solved numerically by a stable and convergent semi-implicit difference method, 
while the counter-propagating spatial evolution is solved by implementing the shoot- 
ing method. We find an enhanced characteristic time scale for idler emission in 
the second-order correlation functions, consistent with the superradiance timescales 
predicted by the analytical method in Chapter 3, and observed experimentally. 

In Chapter 5, the correlated two- photon state derived in Chapter 3 is used to 
investigate the spectral effects on DLCZ protocols involving entanglement generation, 
swapping, and quantum teleportation. We analyze the performance of the protocol 
using, photon-number resolving and non-resolving photon detectors. We find that 
a more genuine and high fidelity protocol requires a source with reduced frequency 
space entanglement. 
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In Chapter 6, we present the analytical results on the efficiency of light-frequency 
conversion in a diamond atomic configuration. We find the optimum efficiency as a 
function of optical density. We find the maximum conversion efficiency by studying 
parametric coupling windows that are created by strong pump fields, and provide 
numerical solutions for the pulse conversion. 
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APPENDIX A 



DERIVATION OF A SCHRODINGER WAVE 
EQUATION FOR SPONTANEOUS EMISSIONS 
FROM A CASCADE TYPE ATOMIC 

ENSEMBLE 

In this appendix, we derive the Hamiltonian for the cascade emission (signal-idler) 
from a four-level atomic ensemble. We use Schrodingcr's equation to study the 
correlated two-photon state from a two-photon laser excitation. Apart from the ro- 
tating wave approximation, non-rotating wave probability amplitudes are introduced 
to take into account the proper frequency shift. The adiabatic approximation on 
laser-excited states is used to simplify the atomic dynamics and solve for the signal- 
idler probability amplitude. 

A . 1 Hamiltonian and Equation of Motion 

Consider an ensemble of N four-level atoms interacting with two classical fields and 
spontaneously emitted signal and idler photons as shown in Figure 3.3. These iden- 
tical atoms distribute randomly with a uniform density. Use dipole approximation 
of light-matter interactions, —d-E where E is classical or quantum electric field, and 
include non-rotating wave approximation (RWA) terms in the interaction of quantum 
fields, the Hamiltonian in interaction picture is 
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Vj{t) = 

N N , N 

/i=i ij,=i /i=i 

N 

-rh.c] -ihY,9k. [ek.^M^.e~'-'''^'''■'' - ^^^^^^^ -4^ 



N 

^^'.E03).(O|e^^^* + |O)^(3|e— 3*], (A.l) 



where the time dependence of laser frequency is absorbed into interaction terms of 
signal and idler fields. Single photon detuning Ai — oua — oui, two-photon detuning 
^2 = + — ^2, and = UJ2 ~ W3. Rabl frequencies are Qa = m\d\\0)S{ka)/h, 
flh = {2\\d\\l)S{kb)/h, and coupling coefficients are gks = {S\\d\\2)S{ks)/h, g^i = 
{Q\\d\\?>)E{ki)/h. The double matrix element of the dipole moment is independent of 



the hyperfine structure, and £{k) = yi^^- Polarizations of signal and idler fields 
a-re ^kiXi ^^^^ direction of dipole operators are dg-, di. 

In the hmit of large detuned and weak driving fields, Ai » that is discussed 

in Chapter 3.2, we consider only single excitations and ignore the spontaneous decay 
during the excitation process. The state function can be written as 
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\m) = 

N N N 

£{t)\0,yac) + vac) + ^ S^(i)|2^, vac) + ^sim^^, 

/U=l /i=l /1=1 fes,As 

W iV 
ks,Xs /' = ! 'I'l -A, /' = ! 



AT AT 



+ E E ^n^)|2.' + E E E ^^7W|3., 3., 1,^,,^, 1,- ,,^) (A.2) 



" V 



where |vac) is the photon vacuum state, s = {kg, Xs), i = {ki, Aj), |m^) = \iTifj)\0)f^~^, 
m = 1,2,3 and |3^,3^) = |3^)|3j^)|0)f^~J. The probabihty amphtudes coupled 
from rotating wave terms in the Hamiltonian are S{t), Afj,{t), B^j_{t), C^{t), Ds^i{t), 
which indicate the complete cycle of single excitation process from the ground state, 
intermediate, upper excited state, intermediate excited state with emission of a signal 
photon, and the ground state with the signal-idler emission. Note that the states 
underlined are coupled through non-RWA terms that describe a transition from upper 
excited state to intermediate one by absorbing a photon for B^{t) and C^{t) and a 
transition from the ground state to the intermediate one by emitting a photon for 
Cf{t) and E{t) or Ci^{t) and C^^O- Apply the Schrodinger equation ih§-^\'il){t)) = 
Vi{t)\'ilj{t)) , and we have the coupled equations of motion. 
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iA„ -- 



^ " 2 2 



HI 



s 

i 



— % 



v<\x 



v>\x 



(A.3) 



The Lamb shift for the atomic transition |3) — > |0) with the optical frequency Us 
and spontaneous decay rate F is du-^[P.'V.{u — cus)"^— P.V.(c<; + cus)"^] that can 
be identified partly within the substitution of these non-RWA terms. We substitute 
Cf into £, C'^'l into C^, and Bi^ into C^, and they are 
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,'0* . . - . . /■* 

X [7r5(cc;i + CU3) - iP.V.(cc;i + ^3)"^] 

= ^ ^ e-'^'^-^'^A^ + z7V£ duji^-PN.{uji + 6^3)"', (A.4) 



i 

- \9i?\{^i ■ d*)? t dt'e'^'''+'^'^^''-'^S{t')\^ l^sit') + e'^'-^^^'^-^^^Cn^') 

i 

du,^VN.{ui + u^r' - Y \9i\'m ■ d*)\'x 
f (^i'e^('^'+'"3)^*'-*) Ve^^*'(""*''-""*'')c;(i'), (A.5) 

Jo _L 

where we have used the symmetric property of ^tixiC) — ^^^(0 = Z^ibiPK^i ' 
'^DP /o* (it'e*('^^+'^3)(«'-t)g«fci-{jv-?>) [21]. The spontaneous decay rate for the idler tran- 
sition is Fj = g^'l^^ja [116, 30], and the same thing for signal transition Vg = g^^j^^l 
that 

= iC^ j dws^P.y.{ujs + UJ23 + (A.6) 

It is now clear the contribution from non-RWA terms to the Lamb shift of the 
idler transition resides in S and C^, which are proportional to N and N — 1. The 
difference of the level shifts then gives rise to — Jq°° duj -^F.V .{u +ujs)~^ , and the other 
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part can be derived from substitutions of RWA terms. The signal transition has the 
same effect as shown in C^. The frequency shift due to dipole-dipole interaction also 
appeared in that has the contribution of interactions from other atoms. There 
also will be contributions from RWA terms, and we will show the complete expression 
for collective decay rate and frequency shift. 

Define Cs,qi — Yl^i.^s^~^^' ^''i substitute Ds,i into C^, and we have 

+ i{N - l)Cs,q^ / duj^-^VN.{ui + a;3)"' 
[p.V.(a;, - oj^)-^ + RV.(a;i + cus)"^] [ J] e'^^^-^'>''^Cs,H - Cs,qi 

POO p 

+ i{N - l)Cs,q, / dj^~Y>N.{uJi + a;3)-^ (A.7) 

Jo 

Renormahze the Lamb shift (last two lines in the above) and use J^^. Cs^^.e'^''^" = 
NC^ then we have 



[ dt'e- 
Jo 



'2^ 
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+ V 



X 



fOO -p 

X / duji — 
'27r 



1^ 



Stt 



dVti[l - {ki ■ di) 



21 Ts 



1^^^ g«(9'i-fci)-*V _ gi(g'i-fei)-r> 



+ (p d^ill - {k ■ dif 
\ki\=k3 on J 



^i{ki-qi)-r, 



E 

9^ 



(A.8) 



Due to the summation of exponential factors from the above, the couphng from 
the other modes (^^ is significant only when q'i — h — qi, so finally we have 



+ iSuoiCs,q, (A.9) 
where the collective decay rate is [22] 

^{Nfl + 1) ^ y 1^ / dn^[l - ih ■ Inj^Y ^^'-'^■''^-'^ . (A.10) 
and the collective frequency shift expressed in terms of the continuous integral over 

a frequency space is 
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r. r 1 



27r 

oo -p 

duji — 
27r 



(A.ll) 



The geometrical constant for a cylindrical ensemble (of height h and radius a) 



IS 



6(Ar- 1) 



where if = k^h and A = k^a are dimensionless length scales, and circular polarizations 
are considered [22]. Ji is the Bessel function of the first kind. 

The alternative way to express the collective decay rate and shift is [21] 



r r r 1 

^ = ^(7V^ + 1) ^ ^- J] F,.(A;3r,.)e-^H^^.-.), 

where 

= ^{[1 - (P • r^,rf^ + [1 - 3(p • f.,)^](^ - ^)}, 
G.,(0 = - ip ■ r.,rf-^ + [1 - 3(p . f.,)^](^ + ^)}, 

and note that ^ = k^Vap- 

A. 2 Adiabatic Approximation 

Under the conditions of large detuned laser excitations, we may use the adiabatic 
approximation to eliminate the laser-excited states and solve for the signal-idler prob- 
ability amplitude. Before proceeding to the adiabatic approximation, we solve C^^g. 
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(A.13) 



(A.14) 



(A.15) 



first and substitute it to solve S^. 



Let Bk^+k^ = ^^e-^Cka+koy^-.B^ and A^, = E^e-^'^^^'^A^, we have 



ks,\. 

r 



s 

X 

27r 



P.V.(a;,-a;23- (A.17) 

where the Weisskopf-Wigner approach is used to derive the decay rate for the signal 
transition, and in conjunction with the result of C*^, the Lamb shift is also derived as 
the difference of level shifts that da;s|^[P.V.(a;s — U22, — ^2)"^— P.V.(u;s + ijJ2z + 
^2)"^]. We then renormalize it and apply the adiabatic approximation. 
When the detunings are large enough that 

|A,|,|A.|»^,^,^. 

We can solve the coupled equations of motion by adiabatically eliminating the 
intermediate and upper excited states in the excitation process. The adiabatic ap- 
proximation requires that the driving pulses are smoothly turned on, and we will 
show under what condition of the pulses that the approximation is valid. 

First we use integration by parts to solve the probability amplitudes in the adi- 
abatic approximation (zeroth order) and their first-order correction. Note that we 
allow time- varying Rabi frequencies. 



Ill 



f 



2Ai 



2Ai 



+ 



oo 



i d 



t 



2AI dt 



(A.18) 



(A.19) 



2(A2 + ir2/2) 2(A2 + ir2/2)2 
where higher order terms involving a second derivative of the fields are neglected due 
to their feature of slow variation. The initial conditions are used in the below, 

d 



B,,+,,(-oo) = Afc„(-oo) = 0, —(na{t')£{t') 



0, 



d 

dt' 



d 



-{^l,{t')AkSt') 



= 0. 



With conditions in the following (i) to (iii), 



(i) 



(ii) 



(in) 



dt 



«1, 



dt 



fJn,{t)A,St) 



«1, 



{A2 + ir2/2)n,{t)Ak^{t) 



we can derive Ak^{t), and £{t) in the adiabatic approximation, 

jVQ„(t)g(t) 



m 

Bka+h{t) 



2Ai 



4Ai(A2+ir2/2) 



1 - 



e-^/-oolf^a(*')Prft' 



2Ai 

iN 



4Ai J_, 



4Ai(A2+ir2/2) 



1 - 



4A1A2 



(A.20) 



(A.21) 



(A.22) 



(A.23) 

(A.24) 
(A.25) 
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where the probabihty amphtude of the first excited state follows the first laser field, 
and the upper excited state follows the products of two laser fields. The AC Stark 
shift is present in the ground state that can be ignored if Ai > iV \na{t')\'^dt' /A. 
This condition is also required for the assumption of single excitations states we 
consider. 

Finally, we have the probability amplitudes associated with the signal Cs,ki {t) and 
signal-idler photons Ds,i{t), 

Cs,ki {t) 

pt piV ^ ^ ^ 



,^(_£|_+i5^.)(t'_t") 



^i{ul,-UJi)t' ^i{uJs-U23-^2)t" ^(^^lly (A. 27) 

Note that e~**^^''"'"*^''^'''^-B^(t') does not depend on the atomic index jjL under the 
adiabatic approximation, and A/c — ka + kb — ks — ki is the phase mismatch. 

The above expressions are the main results of this Appendix and we proceed to 
investigate their properties when Gaussian pump pulses are used in Chapter 3. 
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APPENDIX B 

DERIVATION OF A C-NUMBER LANGEVIN 
EQUATION FOR THE CASCADE EMISSION 



In this appendix, we show the details in the derivations of c-number Langevin equa- 
tions that are the foundation for numerical approaches of the cascade emission in 
Chapter 4. First we describe how to quantize the free electromagnetic field [29], 
and we formulate the Fokker-Planck equation for our system using the positive P- 
representation. We derive the Fokker-Planck equations by characteristic functions 
[27], and the corresponding c-number Langevin equations are derived. The noise 
correlations are found from the diffusion coefficients in Fokker-Planck equations. 

B.l Quantized Electromagnetic Field 

To describe the propagating quantum fields in one dimension, we take the approach 
of the reference [29]. Before proceeding, we specify the positive frequency of a free 
propagating field operator in the discrete space. 



where ek^x and ex{k) specify polarizations of the field, and the interchange of discrete 
and continuous space has relation. 




(B.l) 




where creation and annihilation operators satisfy commutation relations. 




(B.2) 
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For the purpose of describing one-dimensional propagating field (paraxial approxi- 
mation), we discretize the space along the propagation (z) and denote r as the vectors 
on the cross section. We then have 

EHz, r)=^f2 E E (B.3) 



n=-M '2^0 ^ X 

7 27rn , , ,^ ,^ 

A;„ = , ujs,n ^ + Kc , Us^kgC, n^-M,...,M 

where L is the length of propagation that is equally split into 2M -|- 1 elements, 
and the center of the interval is 2; = 2;^ = 2M+1 ^^^^ ~ — M, ...,M. kg is the 
central longitudinal mode of the field. Note that the polarization Cfc a with paraxial 
approximation has k ^ kn- 

The next step is to characterize the transverse mode of propagating field, and we 
introduce a set of orthonormal transverse mode functions {fi,k„±) that J2k„± fi,kn±-fj^^n± 
5ij. A longitudinal annihilation operator is defined as 

Cn,i,\ = flkr, I ^k,X, 

which also satisfies commutation relations [cn,i,x, ^1^, j — Snn'hj^xx'- We can sub- 
stitute ak,x = Z)i Cn,i,xfi,k„jL that 

M 



E^z, r1 = ^ E V ^'^''^^'"'^7V ^ ^ ^'^^■'^^r.,xkK.Cn,,x. (B.4) 

n=—M X kn± 

Let the spatial transverse mode function 

Ui{r) = -^^^'''^^'^kKi. , where J d\dzu*Ui = 1 



and we have 

n=-M ^ ° i,X 

An approximation of a single transverse mode can be applied if only single mode 
is collected for the experiment, and a fiat transverse mode can be assumed {ui = -L) 
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if the collected mode has a narrower spatial bandwidth than the mode function. 
Finally, we have 

E^M^ E \/ff^'^'^^'"^^E^"'^«-Aan,A. (B.6) 

n=-M V ^'0*^ A 

For a demonstration of deriving an interaction Hamiltonian and Maxwell-Bloch 
equations, we use a two-state system (|0) and |1)), and the polarization is not con- 
cerned here. The free field and interaction Hamiltonian (interacting with atomic 
ensemble with atoms) is 

i/o = E^-'"^^^Cn + ^l|l)(l|, (B.7) 



n 



N 



M 

-hg J2 E a'^''"^c„e^('=^+'=")"'" + h.c. 

II n,m=—M 



(B.8) 



RWA 



9 



d^ = a^ + a^'\ E = E+ + E- (B.9) 



hy 2eoV' 

where sums over 2M+i s-toms in the cross sections and the index m on raising and 
lowering atomic operators a characterizes the position of the atoms. The rotating 
wave approximation (RWA) is made in the interaction Hamiltonian and slowly varying 
coupling constant g is taken out of the discrete mode sum and is assigned a central 
frequency, which is the narrow band assumption for the field. 
Now we introduce a new operator 

M 



hi^ , ^ = V Cnc'''-" , l = -M,...,M, (B.IO) 

V n=-M 



which satisfies commutation relations — ^ii'^ ^-iid the Hamiltonian can be re- 

expressed as 

I w 

V ^-hgY^ V2M + W'Hk^^^'' + h.c. (B.12) 
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The Heisenberg equation of slowly varying field operators {bi = bie^'*) is 



n 



(B.13) 

and we may use the limit of M — > oo that 

Zm = <2M + l ^ ^ ' ^^2MTi hi ^ Et{z, t) , -i^a;„4V2M + l ^ -c—E+{z, t) 

(B.14) 

where the derivative can be shown from 

^ 2M + 1 2M + 1^^ dzi ^ 

I' V n I' n '■ 

--<^y^Sii'bi' = -c^bi- (B.15) 
dzi dzi 
v ' ' 

In the end, we have 

^^^'^^^^^^^^'^^ = + (B.16) 

and use the limit, 

r 2M + 1 

limM-).oo — ^ % 5(2; - 2; ), 

then we have (define slowly varying atomic operators o"'^'' — (jMi'g-^'ss^^z+i'^stj 



AT TV N: 



L 



N 



zi->z 



//=1 //=1 /i 

The field propagation equation in Maxwell-Bloch equations becomes 

U + Cg-)Etiz, t) = ig* E ^"K^. - -i9*^Y, -o"- (B.18) 

B.2 Positive P -representation 

The phase space methods [32] that mainly include P-, Q-, and Wigner (W) rep- 
resentations are techniques of using classical analogues to study quantum systems. 
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especially harmonic oscillators. The eigenstate of harmonic oscillator is a coher- 
ent state that provides the basis expansion to construct various representations. P 
and Q-representation are associated respectively with evaluations of normal and anti- 
normal order correlations of creation and destruction operators. W-representation 
is invented for the purpose of describing symmetrically ordered creation and destruc- 
tion operators. Since P-representation describes normally ordered quantities that 
are relevant in experiments, we are interested in investigating one class of generalized 
P-representations, the positive P-representation that has semi-definite property in the 
diffusion process, which is important in describing quantum noise systems. 

Postive-P representation [35, 94] is an extension to Glauber-Sudarshan P-representation 
that uses coherent state (|q;)) as a basis expansion of density operator p. In terms 
of diagonal coherent states with a quasi-probability distribution, P(q;,q;*), a density 
operator in P-representation is 

p= I \a){a\P{a,a*)(fa, (B.19) 

JD 

where D represents the integration domain. The normalization condition of p, which 
is Tr{p}= 1, indicates the normahzation for P as well, P{oi, a*)(Pa — 1. 

Positive P-representation uses a non-diagonal coherent state expansion and the 
density operator can be expressed as 

p= /" A(a,/5)P(a,/3)ci/i(a,/3), (B.20) 

JD 

where 

dp{a, /3) = (fad?^ and A(q;, /3) = , (B.21) 

and in non-diagonal projection operators, A(a,/3), makes sure of the normal- 

ization condition in distribution function, P{a,(3). 

Any normally ordered observable can be deduced from the distribution function 
P{a, 13) that 
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((a^)'"a") = / ^VP{a,l3)dii{a,f3). (B.22) 
Jd 

A characteristic function Xpi^a,^^) (Fourier-transformed distribution function in 

Glaubcr-Sudarslian P-rcprcscntation but now is extended into a larger dimension) 
can help formulate distribution function, which is 

XpiK,X^)= f e^^«"+^^^^P(a,/3)ci/i(a,/3). (B.23) 
Jd 

It is calculated from a normally ordered exponential operator E{X), 

Xp(Aa, A^) = Tr{p^(A)}, ^(A) = e'^^'^' e'^"\ (B.24) 

Then a Fokker-Planck equation can be derived from the time derivative of char- 
acteristic function, 

by Liouville equations, 

| = ^I""'l- (B.26) 
In laser theory [27], a P-representation method is extended to describe atomic and 
atom-field interaction systems. When a large number of atoms is considered, which 
is indeed the case of the actual laser, a macroscopic variable can be defined. Then 
a generalized Fokker-Planck equation can be derived from characteristic functions 
by neglecting higher order terms that are proportional to the inverse of number of 
atoms. It is the similar to our case when we solve light-matter interactions in an 
atomic ensemble that the large number cuts off the higher order terms in characteristic 
functions, which we will demonstrate in the next subsection. 
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B.2.1 Hamiltonian 

The Hamiltonian is in Schrodinger picture, and we separate it into two parts where 
Hq is the free Hamiltonian of the atomic ensemble and one dimensional counter- 
propagating signal and idler fields, and Hj is the interaction Hamiltonian of atoms 
interacting with two classical fields and two quantum fields (signal and idler) . Dipole 

— * — * 

approximation of —d • E and rotating wave approximation (RWA) have been made 
to these interactions. Similar to the previous Appendix, we have 



H — Hq + Hi , 

3 M M 



i=l l=-M 
M 



l=-M 



1,1' 



+ fujiY^ dlidi,i + h^ujii>dlidi^i> , 



l=-M 
M 



1,1' 



l=-M 

M 

l=-M 



(B.27) 
(B.28) 



(B.29) 



where a^,^ ^ = E^^ \m),{n\ _ , ^ /„(i)dio^(/c„)/(2;i), and is 



rij,=zi 



slow varying temporal profile without spatial dependence (ensemble scale much less 



than pulse length). Qs = d23£{ks)/h, £{k) = ^Jhw/2eQV and 



mL 
2M+1' 



— M, M, and L is the length of propagation. Note that the Rabi frequency is half 
of the standard definition. 

The normally ordered exponential operator is chosen as 
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I 

E^(^X) — e*'^19^01 g'-^18^12 g*-^17^02 g^-^ie^lS g*''^'l6^03 g*-^14^32 g^-^lS^ll g*'^12^22 g^-^ll^lsg'-^lO^L 

g^4'5'03e*-*'8°'l3g«^7°"02g'-^6'5'l2g'^5'^W (B.30) 

Aside from the atom-field interaction |f = ^ [H, p] , when dissipation from vacuum 
is considered (single atomic decay) , we can express them in terms of a Lindblad form 
where we have for the four-level atomic system, 

l=-M n 

^ 2 1-12 P'-'l2 '^12 '^12 P P^12 ^12 J 

+ ^f[^^'JpKt' - K^^<ip - p^'.^'^<i\ 

The characteristic functions can be calculated, 



X = Tr{£;(A)p}, (B.32) 

+ (|).-.+ (|),. (B-33) 

(|).-. = TV{B(A)^4|ff.-.,pl}, (|),, = (B.34) 

where i?o = Ha + -f^L, -f^A is the atomic free evolution Hamiltonian, Hj^ is the 
Hamiltonian for laser fields, and Ha_l = Hj. Now we continue to derive the time 
derivative in each part of characteristic functions. 
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B.2.2 Chciracteristic function - atomic pcirt 

The atomic part in characteristic function is deduced from (|^) 



(|), = TVmA)(|),}, 



I 

so various components in (^)^ are 



Tr{E{X) J2 ^iiP} = E T^im^nP} 

I I 



^ ^ a(«Ai9) 0(«Ai8) 5(^Ai6) <9(2Ai3) 

"^{EiX) E ^22p} = Et^'^eT^TTT-^ + ^Ar + iAio + . J ;X, 

^ ^ a(2A6) <9(^A7) d(i\io) d{t\i2) 

d d d d 

Tr{E{X) J2 P-L} - E[^^-a(^ + '^^^0^) + + mu)^'''' 

^ ^ 

E ^^33} = Et^^iS ;.... x +^Ai5 . ^ -^Al4 . . + , A ix 
^ ^ d[%\x%) d{t\i5) d{iXu) d{tXn) 

(B.36) 

where the subscript I on the bracket reminds us the derivatives inside the bracket 
operate on Zth component of the characteristic functions. 

B.2.3 Chciracteristic function - field part 

The field part in characteristic function is deduced from 
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(|), = TVmA)(|)J, 

dp X 

(^)l ^ ^[^l,P\^^\- i^s{al^ias,iP - pa^iKi) - i^i{a{iai,ip - pa{iai,i) + 

I 

[ - iuvi{al^ias,vp - pa\ias,v) - iuii,{a\iai,i,p - pdlidi^i') , (B.37) 

1,1' 

and various components in (^)^ are 



^2 d ^ 



and 



Tr{E{\)J2ali<iP} 
I 


I 


Tr {E{\)J2pal,ids,i} 
I 




Tr{E{\)Y,^^,A^,lP} 
I 


= Ei 


TT{E{X)J2pdld.,i} 
I 


I 



92 5 1 



(9^ (9 

^-^27^7TVT]a, (B.3^ 
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Tr{^(A) j:u,,dld,M = E-'''[a(iAi)a(iAD+^^^'a(iAf)l^' 



(9^ (9 
Tr{E{X) Y,^,ipdlds,i'} = E^'''[7)77):Tw7in + ^^4 7^]X, 



B.2.4 Characteristic function - atom-field part 

The atom-field interaction part in characteristic function is deduced from 
and we denote part (a) for the classical field interaction. 
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M 



■dt' 



ih 



l=-M 



(B.40) 



and various components in (^)^l^ are 



Tr{E(A) J]e^^-'4\p} = 



d 



d 



d d d 



d 



diiXa) d{iX 



d 



11 



d{iX 



d d 
2^-^5^;:7TT — T + iX^Nz — iX-i- 



12 



'd{iX 



13 



d{iXe) 



r\ r\ O 



d 



d{tXs) 



d{iX 



15j 



d{iXig)U 



X, 



T:t{E{X)Y,e-""'''al,p}^J2''~'''' 



d 



id{iX5 
r d 



X: 



Tr{£;(A) Epe-''"^'<} = IV{£;(A) ^ e''"^'4p}*A^o-A« a*o-a,s.a*«-a,. A-^-A.a, 

; I Ag-<->— Ai5, Aj^<->— Aii,Aj2<->-— Ai2,A5;3-H-— Ai3 

(B.41) 



where a correspondence that we denote as C later is Ag ■<-> — Aig^ Ag ■<-> — Aig^ Ay •<-> 
— Ai7^ Ag -H- — Ai6, Ag -H- — Ai5^ A^i -H- — Aii^ A*2 -h- — A12, A^g -h- — A13, can be observed 
to help calculate the characteristic function. Also 
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I 



9 5 9 9 9 

(iA6)(iAio) , + (iA7)(iA6) , + iAe ^,.. , - iAe ^,.. , + 1X7- 



d{iXio) ^(iAy) 9(iAi3) 9(iAi2) ^(iAg) 



9(iAi2) 9(iAii) 9(iAio) 9(iAi4) 

X, 



+ iAioe^^"-'^" . + ((iAio)(iAi4)e'^"-^^" + e^^^^-'^^^^ ^ 



9(iAi6) 9(iAi8) 



X, 



I I 

Tr{E{X) J2 Pe""'''^[2} = T^{E(X) J] e-'^^^'af, p}^, (B.43) 
I I 
and the atom-field interaction characteristic function for quantum fields, which we 

denote as part (b), is 

«=-M 

+ giV2M + lagai,ie''^'^ + h.c] , p] . (B.44) 
For the part of fields only, 



TT{E{X)a,,,p} = [q^^]x, Tr{E(A)pa,,,} = + tX[]x, 

Tr{EiX)alp} = + iX^jx, Tr{E{X)pal} = [^]x, 

Tr{E{X)ai,ip} = Tr{E{X)pai,i} = [—^^+iX^^]x, 

Tr{E{X)dlp} = + iX[]x, T:^{E{X)pal} = (B.45) 
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and for the part of atomic operators associated with signal field, 



Tr{E{X)alp} = 



. d 



d{iX,) 
d 



{iX 



d 



10 1 



iXi 



d{iXr^ 
d 



+ 



(9(zAi2^ 
d 



iAi2— iAii 



d{iXio) 
d 



d{iX 



X, 



14;-!' 



Tr{^(A).^,p} ^...^^^^^^ , ^^^^^^^ 

IV{^(A)p4} = Tr{i?(A)42P}o 
T,{E{X)pdi^} = IY{£;(A)4p}^, 



+ «Ag 



— 1 



X, 



(B.46) 



and for the part of atomic operators associated with idler field, 



Tv{E{X)~a%p} 

{iXe){iX8f^7T^ + (iA5)(iA6)(iA8)^7TVT - (iA5)(iA8)( , -iXg, 



d{iXs) ^(iAe) d(iXi3) d{iXg) 



d{iXio) d{iXu) d{iXio) ^(iAie) 

(^A5)(^Ai4)e^^"-^^"^-^ - (^A5)(^A9)^7^ - (^A7)(^A8)7^ - iAy- ^ 



^(iAis) ^(iAs) ^(iAe) 9(iAio) 

<^ f) 

- {iXr){iX,)—— + iX,N, - (iAg)^-— + iA8[-(iA9)— - +iAi6e^"" 



d{iXn) d{iXi9) d{iXi5) d{iXi3) d{iX 

^ +2^Alo^7:^-2-^]+^AMe''"^7^1 X, 



w 



d{iXi2) d(iXio) d(iXn) d(iXn) 

d 

■ X, 



Tr{i^(A).^3P} 
IY{£;(A)p4} = Tr{E{X)al,,prc, 

Tr{E{X)pal,} = Tr{^(A)4Uc- (B-47) 
B.2.5 Characteristic function - dissipation part 

We calculate the characteristic function from (§7)^^ up to the second order of various 
A's (where we denote (2)) that account for drift and diffusion terms in Fokker-Planck 
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equation. Below we drop the summation over spatial slices I, which we will retrieve 
later, 

= Tr{[(iA8)(iAio)<Ti2 + (iA6)(iAi4)(Ti3 + (iA8)(iAi5)aio+ 

(iA8)(iAi6)aii + (iA6)(iAi7)<Tio + e-'^'^an - iXige-'^^'aw 

+ {iXe){iXis)au - iAge-^^^aia - iXee-'^''ai2]E{X)p}^^^ (B.48) 

where various properties of tracing can be found in previous sections, and the one we 
did not have before is (up to first order) 



Tr{ai3^(A)p}(^) 

d d 

= [^Ai5 - iXi6{——— - ) - iXi8- 



d{iXig) d{iXii) d{iXi3) d{iX 



14 



d{iXs) oiiXe) 

Put everything together, and for the dissipation of first laser transition we have 



1 1 

7oiTr{E(A)[(Joip(jJi - -anp- ^P^n]}^^^ 



7oii 



iAs d 



iX 



19 



d 



iXa d 



iAis d 



iXii d 



2 d{iX^) 2 a(iAi9) 2 ^(iAe) 2 d{iXxs) 2 d{iX^ 

- «Ai37^7T^; ^ + {lXi3){lXi8) ^, . + (^Al3)(^Al6J 



2 diiXie] 



d{tX 



13 



d{iX 



18 



diiXe 



d d d d 

+ (iAi3)(iAi6) . + (iAi3)(iA8) . + (iAg) (i Ais) ^ttt-^ + (iA6)(iAi6) 



d{iX 



16j 



d{iX8) 



d{iX 



14 



+ (^A8)(^Al6) . , + (^A6)(^Al8) , , H r — J x- 

9(zAii) 9(zAi2) 2 9(zAi3) 

And for the second laser. 



d{iXw) 
(B.50) 



Tr{E(A)ai2paL}(') = Tr{ [{iX^^){iX^^)a2o + {iX^&){iXi^)(T2i + e^^^-^^^^ (^22 - iX^^a^o 
- iAi4a-23 - iAi8a-2i + (iAi9)(iAi8)(72o ) -E'(A)p}^^^ (B.51) 
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The above requires 



Then we have 



7l2Tr{£;(A)[(Ti2p(Tl2 - ^<3-22p- = 

,. iAe d iAi8 d i\r d iXu d iXw d 
7i2[- 



2 diiXe) 2 d{iXis) 2 ^(^7) 2 a(iAi7) 2 d{iX 



^^T^x^ + (*Ai3 - ^Ai2) . . + (^A5)(^Al9)— — — H x- 

2 d{iXu) d{iXi2) o{iXi2) 2 a(zAi2) 

(B.52) 

And the dissipation for the signal transition, 

Tr{i?(A)a32p4}^'^ 

Tr{[(iA8)(iAi6)o-22 + (iA9)(iAi5)5-22 + (iAi4)(iAi5)5-2o + (iAi4)(iAi6)o-2i 
+ e^^""^^i2(-iAi70-2o - iAi46-23 + 0-22 - iAi80-2i + {iXig){iXi8)^2o)]E{X)p}, 

(B.53) 

so we have 



732Tr{£;(A)[0-32pCTj2 - ^^22P - ^P^22]}^^^ = 

,. iAe d iXi8 d iXr d iXn d iXio d 

732 [ 



2 d{iX(i) 2 9(iAi8) 2 a(iA7) 2 9(iAi7) 2 d{iXw) 

^Al4 d ..X <9 /.xN/.xN d /.^x/.^^ ^ 

+ («Aii - 1X12) ^,.. . + (?A5)(^Ai9) + (^A8)(«Ai6)- 



2 9(iAi4) ^ ^a(iAi2) ^ -^^^ ^^^9(^12) ^ ^"^5(iAi2) 

And for idler transition, 
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Tr{[(iAio)(iAi8)a3i + {i\w){i\i4)^33 + {i>^w){i>^n)^30 - i><ioe~''^^^a32 
e-'^"(a33 - iAi6(T3i + (iAi6)(zAi9)a3o - iXi5aso]E{X)p}^^\ (B.55) 

The above needs 



then we have 



7o3Tr{^(A)[(Jo3po-j3 - ^0-33^ - ^^0-33]}^^^ 

iXs d iXiQ d iXg d iXirj d iXiQ d 

= 703[ 



2 d{iXs) 2 d{iXi6) 2 d{iXg) 2 ^(iAis) 2 ^(iAio) 

5.5 Stochastic Differential Equation 

A distribution function can be found by Fourier transforming the characteristic func- 
tions, 



/(«) = / - / e-'""\(A)dAi...dA„, (B.57) 

then 



9f 1 f f ^-^s■^9x 



, , , ... , e-"^-^^dX^...dXn. (B.58) 
dt (27r)" 7 7 9t ^ ^ 

If ^ = iX^ ) ' integration by parts and neglect the boundary terms, we 
have ^ = —-g^^a^f where a minus sign is from iA/3. Correspondingly, if ^ = e*^^, 
we have ^ = e ^'"''^ . 
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B.3.1 Fokker-Planck equation 



Let 

1,1' 

then we have for the atomic part, 



-ia„[^(-<.^) + ^(-„i) + g|-(-„;„)l 

where C is ctg •<->■ ctig, ctg •<->■ ctis, cty •<->■ ctiy, ctg •<->■ ctie, ctg •<->■ q;i5, q;*q •<->■ q;i4, aij^ •<->■ an, 
q;^2 o;i2, Oils •<->■ Q;i3, 0:2, ctg •<->■ a^, and c.c. is complex conjugation. Also for 

the field part, 



_d_ I 
d I' 



d 
d 



+ iuji 



d 



d 



dai 



Ta^ + iujii> 



_d_ 
dai 



_d_ 



(B.61) 
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The atom-field interaction part (a) is 



A") _ 

'-'A-L — 



+ e"4^^4^)«)] - ^fifee-^*^''^'-^'^^*[al8 + -^{-a[,)] + (c.c. with C"), (B.62) 

and let = ^a'-l s + ^a-l/' which are the terms for signal (S) and idler (I) 

parts, 



,(6) 



V2M + le-''''^' I 



d 

dai 



:-«9) - 



da^Qda^Q 



(«lo) 



+ -^(-< + Q;l2) + e ^ ^(al4)]4 + i^:V2M + le^'=»^' 



+ (c.c. with C) 



(B.63) 



and 
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a3 a3 op. 



tts) + TT^e ^"13 (-a 



-2ay + ^e"4T(-al,)]al 

+ i^*V2M + le-**=^^'(a^)(4 - ^) + (c.c. with C). (B.64) 

The dissipation part Cgp can be derived accordingly and the above equation, which 
involves higher order derivatives (third order and higher), is neglected. The validity 
of truncation to second order is due to the expansion in the small parameter l/N^. 

If the Fokker- Planck equation is 

where A and D are drift and diffusion terms then we have a corresponding classical 
Langevin equation 



^ = A^ + T^,^=A, + T, (B.66) 
with a correlation function (Far^) = S{t — t')Da/3. So we have according to various 
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- igiV2M + le"'^''a[ea[ + Vl 



14"3 



- i^:V2M + le^'=»"'4o«i + r'ls, 
a'li = -7o3an + 732al2 - igsV2M + le-'^'^' a{^a\ + ic/,V2M + le^'^^^'^oO;! 
+ i^y2M + le^*^^^'454 - i^*V2M + le-^*^*^'44 + Tn, 



2 



2 

+ i^;V2M + le^'=^^'(42 - 4i)4 + i^iV2M + le*'='^'4X + 
I' 

a[ = -iujia[ - i ^ uu^a^ + ig*V2M + le-'^^"a^g + r[, (B.67) 
I' 

where 72 = 712 + 732- We postpone the derivations of diffusion coefficients after the 
scahng is made in the next section, and note that the complete equations of motion 
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are found by making complex conjugate of the above with correspondence C and 
changing Langevin noises correspondingly, say Fg Fig. 

B.3.2 Slowly varying envelopes and scaled equations of motion 

Here we introduce the slowly varying envelopes and define our cross-grained collective 
atomic and field observables, then finally transform the equations in a dimensionless 
form for later numerical simulations. We note that 



i^^^ii'Oi'i = ~ i^^^u'O^i = "'^^'^1' (B.68) 

and Uq — Nz — ct^g — ~ ^Ai- Define slow varying observables that 



1 1 

^, (~ +\ — ^-ikaZi+ikbZi+ioJbt+ioJat ~ +\ — ^-ioJat+ioJst+ikaZi-ikiZi 



~ 1 ~ 1 

ai2{z,t) = —a[r^,ais{z,t) = -^a[^, 

1 



au{z,t) = —a[^e-'^'^'''+^-'^'e''''^''-''"'''-'''^'' (B.69) 

^z 

where e*^^^ = e^kazi-iktzi-ikizi+ik^zi _ ^^g^ ^j^^ g^j^^ variables, 



E;{z,t) = -f^V2M + la[e-"^'\ E^{z,t) = -p-V2M + la[e''''', (B.70) 
di/n di/n 

where we use the idler dipole moment in signal field scaling for the purpose of scale- 
free atomic equation of motions, so we need to keep in mind that in calculating signal 
intensity or correlation function, an extra factor of {di/dsY needs to be taken care of. 

Wc choose the central frequency of signal and idler as ujs = 0023 + A2,UJi = 
where Ai — oua — oui and A2 — ooa + ouj, — 002- With a scaling of Arecchi-Courtens 
cooperation length [115], we set up the units of field strength, time, and length in the 
following, 
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n ' d,/h ' - T, V 2heo ' 260 T,d,/h- ^""-''^ 

Compared with optical density and superradiant time scale, we have (in terms of 
single atomic decay rate 7) 

Now the slowly varying and dimensionless equations of motion with Langevin 
noises in Ito's form are 



—55 = (iAi - ^)55 + i^a{ao - Sis) + iillar - laiQE^ + J^5, 
—56 = i(A2 - Ai + ^^ )56 - iJl*57 + iJl6(5i3 - 5i2) + iasEte-'^^^ + J'e, 

—57 = (iAs - y )57 - ^f^aSe + i^bf^b + iagE^e''^''^ - iaio-E+ + J^7, 
— q;i3 = -7oiq;i3 + juau + iO,aaig - iQla^ - iQ^aig + i^lae + J"i3, 
— q;i2 = -72Q;i2 + if^fcOiig - iQlae + iauE^ e''^'"'' - iaioE'e'^''' + Tvi, 
—an = -7o3aii + 732ai2 - iau-E^^e"*^'"'' + iaioE" e'^'''' + iai^Ef - iagE~ + J^n, 
-58 = -(iAi + ^°'^^°' )58 - i^lag - in^au + ia^Eje'"^^' + ia^gE+ + Jg, 
-Q!9 = -— Q!9 - iOaCts + iarE'e'' + ^(ao - + -7^9, 



^Qg -|- ^2 „ • A r. I 

au = -(iA2 H )au — i^l^s + ^(0:12 - Oin)E~e'^ ^ + ianE^f + Tu, 



dt 

(B.73) 

and field propagation equations are 

(| + ^)£^^ = ^59 + ^1, (B.74) 
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where is a unit transformation factor from the signal field strength to the idler 
one. For a recognizable format of the above equations used in the text of Chapter 
4, we change the labels in the below, 



as TToi, aQ -H- 7ri2, ar -H- 7ro2, as '^n, ag ttos, aio -H- 7:32, an -H- ttss, 

ai2 -H- 7722, ai3 -H- TTii, 5i4 •H- 7r|2, ai5 Til^, 5i6 •H> Trja, 5i7 •H- 7r^2) 

ai8 ^ 7r|2, ai9 O tt^^, (B.75) 

where tt^ is the stochastic variable that corresponds to the atomic populations of 
state \i) when i = j and to atomic coherence when i ^ j. Note that the associated 
c-number Langevin noises are changed accordingly. 
The Langevin noises are defined as 



1 1 

T {■r — r' i-\ — T^l „-i0Jat+iuJ2,t+ikaZi-ikiZl 

-''Z -''2 

= -jjT\/2M + le-^'"=*ri, = -^V2M + le^'^'*r'i (B.76) 

where other Langevin noises can be found by using the correspondence similar to C", 
for example, •H- J-'ig. 

Before we proceed to formulate the diffusion coefficients, we need to be careful 
about the scaling factor for the transformation to continuous variables when numerical 
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simulation is applied. Take (J-qJ-^) for example, 



{T,{z,t)T,{z\t')) 

_}_'iki,zi+ioJi,t -ikaZii+ioJat' /'pl'pl' 



1 

z 

1 



S{t - t')5{z - z') 



[iinaTc)ae + iaioiE+/E,)] —8{t - t')TJ(z - z')L. 



2M + 1 

L iV. 1 



Lr N 



= ^ [z(0„re)56 + ia,o{E+/E,)] ^S{t - t')T,5{z - z')L, 



(B.77) 



where we have used limM^oo ^^5/;. = 5{z - z'), 2M + 1 ^ and ^ ^ is the 
cooperation number. Then we have the dimensionless form of diffusion coefficients. 



{M^, t)Mz', t')) = j^D,,,6{t - t')S{z - z') (B.78) 
L>6,5 = {i^oPi% + «5io£;+] . (B.79) 



The dimensionless diffusion coefficients are 
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(i) £)5,5 = -i2fia55; D^^Q = iiO,aaQ + aioE^); D^j = -iO,aar; 
-05,8 = ii^a^s + - 5i3)£'+); Ds^g = -ii^aag + 55£;+); 

-^5,11 = -iSie-Bf ; -Ds.is = iotiQEf] Ds^u = -i5i8-B^; .05,19 = 7i2ai2; 

(ii) L'6,6 = -i20b56; £'6,8 = -i^bOis] D^^io = -i^b^io; 

-06,13 = -i0.laY + 7oi56; -D6,i6 = -iot^E' + 7oi5io; i^6,i8 = 7oi5i2; 
(m)L'7,8 = -ia^El; Dr,9 = -iarE:^; 

{iy)D8,9 = -iasE^; Ds,io = ^^^6(5l2 - Sn); i:)8,ii = if^bftM; 

-^8,12 = -i^bOiU] -^8,13 = -?f^aS9 + iSi9E+ + 7oi58; 

-^8,16 = iauE^ - iagE~ + 7oi5ii + 732212; -D8,i8 = iotnE^ + 7oi5i4; 

(v)£'9,9 = -i2agEp, Dg^io = iaioEl; Dg^is = 7325l2; 

(vi) Dio,io = -i2aioE+e~'^'''; Dio,ii = i(nbai6 - oc^Er) + 7o35io; 
-Dio,i3 = -if^feSie; -Dio.M = i^bOiis, - i^lotQ + 703S12; -^10,19 = iaeE^; 

(vii) Dii,ii = iauE+e-'^''' - iawE'e'^^^ + iair,E+ - iagE' + 732512 + 703S11; 
i^ii,i2 = ta,oE:e'^''' - zSi4^+e-^^^^ - 732X12; 

(vm)Di2,i2 = i^^65i8 - iniae - iaioE'e'^''^ + iauE+ e''^''^ + 725i2; 
-Di2,i3 = -ififeSis + i^lae - 7i2ai2; 
(ix)i:)i3,i3 = iQaSiQ - iQ*55 + iQbOiis - i^l^e + 7oi5i3 + 7i25i2; 

I |2 I |2 

(x)D3,8 = f\^taee'^''; ^3,9 = f^^5re^^'=^ (B.80) 

\9i\ \9i\ 

B.3.3 Alternative method to derive diffusion coefficients by Einstein re- 
lations 

Before going further to set up the stochastic differential equation, we show here 
how we derive the diffusion coefficients from the Heisenberg-Langevin approach with 
Einstein relations, and it provides the important check for Fokker-Planck equations. 
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We note here that a symmetric property of the diffusion coefficients is within Fokker- 
Planck equation, whereas the quantum diffusion coefficients in quantum Langevin 
equation do not have symmetric property simply because the quantum operators do 
not necessarily commute with each other. 

The approach involves a quantum-classical correspondence in deriving c-number 
Langevin equations and requires a chosen normal ordering of quantum operators. We 
use the same ordering as we use for deriving Fokker-Planck equations in Eq. (B.30), 

^01> ^12) ^02) ^13> ^03> ^32> ^22, 5"33, 0-32, O-qs, 5"i3, Cro2, Cri2, CTqi, SJ, Og, oj, tti 

and its classical correspondence is 5i9_ai8_...5i. 

We take -Ds.ia = -^13,8 for a demonstration. We first calculate the quantum 
diffusion coefficient, 1*13,8, using Einstein relations where we attach the hat to it, and 
then we can find "^13,8, a classical diffusion coefficient, which is reviewed in Chapter 2. 
Note that in calculating the quantum coefficients, we take advantage of Eq. (B.73) 
where the drift terms are directly corresponded to quantum Langevin equations. For 
clarity, 1^13,8 = ^^11,^13 with 5i3 — > au,as — > d'ls representing a correspondence to 
quantum Langevin equations. The index in the classical variables a represents the 
ordering we choose as defined above, and in various quantum operators a, the index 
represents the atomic levels for atomic coherences or populations. We should find 
"^13,8 = -Di3,8 = -08,13, and the proof is illustrated below by the Einstein's relation, 
Eq. (2.8), 
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V 



13,8 / - 



-7oic^ii + 7i2cr22 + i^a^li " i^lo-Qi - iflba[2 + iflla^ 



;;-t 



(713 



0-11 



d 



+ ^ (anaia) 

= (7oi5"i3) , 



(B.81) 



where the term ^ (5"iia"i3) = ^ (^13) is the drift term of the quantum Langevin 
equation that can be found from Eq. (B.73), 



-ai3 = -(^Ai + 2^1±^)ai3 - ^K^o^ - ^^b42 + ^&l2E;e'^''' + zalE^. (B.82) 
Prom Eq. (2.15), we have 

(^13,8> = 

^13,8^ + {{ -701O-11 + 7125-22 + i^a^li " i^l^oi " i^b^l2 + i^l^i2 
an 



0'13 ) + 



classical counterpart | 



(B.83) 



where classical counterpart represents the last two terms of Eq. (2.15). We have 
used the commutation relations for non-normal correlation functions that 



[5"oi 


5-13] 


= 5-03, 


[5-12 


c^is] 


= 0, 


5"ii, 


^32] 


= 0, 


C^ll, 




= ^01, 



(B.84) 
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and use the correspondence CTqs — > ctg and CTq^ — > ctig in X>i3,8- The rest of the diffusion 
coefficients are confirmed by the method of Einstein relations illustrated above. 

B.3.4 Ito and Stratonovich stochastic differential equations 

The c-number Langevin equations derived from Fokker-Planck equations have a di- 
rect correspondence to Ito-type stochastic differential equations. In stochastic sim- 
ulations, it is important to find the expressions of Langevin noises from diffusion 
coefficients. 

For any symmetric diffusion matrix D{a), it can always be factorized into 

D{a) = B{a)B^{a) (B.85) 

where B BS (an orthogonal matrixiS that SS'^ = I) preserves the diffusion 
matrix so B is not unique. The matrix B is in terms of the Langevin noises where 
^idt — dWl (Wiener process) and {ii{t)^j{t')) — 5ijS{t — t') and the below is just a 
random number in Gaussian distribution with zero mean and unit variance. 

In numerical simulation, we use the semi-implicit algorithm that guarantees the 
stability and convergence in the integration of stochastic differential equations. So a 
transformation from Ito to Stratonovich-type stochastic differential equation is nec- 
essary, 

dxi^ Ai{t,^t)dt + ^Bij{t,^t)dW^ (Ito) (B.86) 

j 

dx\ = [Ai{t, 4) - ^ ^ ^S,fc(i, ^t)-^Bik{t, ^t)]dt 

j k 

+ '^Bij{t,x^t)dW^ (Stratonovich) (B.87) 

j 

where a correction in drift term appears due to the transformation. 

Here we have the full equations with 19 variables in the positive-P representation, 
64 diffusion matrix elements, and 117 noise terms (random number generators). A 
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correction in drift term is underlined and we have (S for Stratonovich) 

= (^) +(iAi - ^)55 + i^a{(io - S13) + ifllar - iaieE^ + T^, (S) 

^5i9 = ( ) +(-iAi - ^)5i9 - in*(5o - Sis) - iVL^fiirj + ia^El + J^g, 

^5i8 = -i{^2 - Ai - i ^°^^^^ )5i8 + if^a5i7 - i^l{ai3 - Sia) 

= (iA2 - 5)57 - i^aae + i^b^s + ia^E^e''^^^ - iaioE+ + JV, 

OT Z 

= (-iAi - '^"^ ^^°^ )58 - ^f^a«9 - if^ftSu + iaQE'e'^^"" + iSigE/" + J'g, 
— 5i6 = (iAi )5i6 + i0.aOti5 + i^^bSio - iaigEfe''^'''' - ia^E' + J^ie, 

OT Z 

— ag = {iE'l) — — Q;g - iOaCts + iar-^s" + ^(ao - ocxi)E^ + J'g, 

C/T ^ — ^^—^ Z 

— ai5 = (^^^^^^72 — 2^*^^^ ^^a"i6 - ianE^ - i{ao - au)E~ + J'ls, 

- ia-jE^ + Jio, 

^5i4 = {-'-E:) +(-zA2 - - i^a^ + ^(Sl2 - 5n)^;e'^^^ 

ar z z 

" V ' 

+ ianEf + ^"14, 
9~ -5701 + 712 ~ . ~ ~ *~ ~ 

T^aiS = 70iQ;13 + 7120:12 + l^^laOH^ - ^iiaO^S - ^i^bttlS 

dr ^ 4 ^ 

+ iQ^Se + J"i3, 

^5i2 = - ? -725i2 + if^bSis - iO^aQ + iau^+e-^^'^^ - iawE^e'^''^ + Tvi, 

^ -3703 + 732 ^ ^ ^ _ ,5^^£;+e-iAfe. ^ i5^^£;-e^Afc. ^ ^5^^^+ 

dr ^ 4 ^ 

- iagE^ + J^ii, 
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d 

—E+ = iag + J^i, 

^E- = -ia,5 + T2. (B.88) 
The Langevin noises are formulated as a non-square form [35, 84] 
Tx^T2^ 0; 
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— 1^(62 - ^63) + \/ ^19,196 + y^^(68 + ^69) + ^^^(CSO + ^Csi) 



-Dig, 16/^ , \ , , \ , /-Dl9,l0, 



+ aF^(^32 + ^^33) + A ^(^34 + ^e35) + A ^(^36 + li.l) 



-Dig, 13. I -t \ I /-^19,ll/(^ , V 
-^^(^38 + <39) + A/ ^^^(^40 + <4l); 



•^6 = a/ - iilh) + \/^C3 + a/ — 1^(^42 + i^43) + \\ ^(^44 + i^A^) 



+ a/ ^(^46 + ^^47) + A/ ^^(^48 + ^^49) + \\ ^(^50 + ^60; 



:^18 = 



18 = y^^(68 - i^g) + ^J^^{U2 - iU3) + V ^18,18^4 + ^J^f^{C52 + iC53 
+ a/ —^7^(^54 + «C55) + a/ —^7^(^56 + ^^57) + A/ —^(^58 + ^69); 



•^7 = a/ -^(^16 - ^67) + a/ -?7^(C60 + ^Cei) + a/ -?T^(C62 + ^^63); 



^17 = a/ —^(^30 - iisi) + a/ —^(^64 + iCdb) + a/ ^ o'^^ (^66 + ^^67); 
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•^8 = ^^(Cl8 - ^69) + ^^(^4 - ^^45) + ^^^(^52 - ^^53) 
+ /^(e72 + ^e73) + + ^e75) + ^^(^76 + ^C77) 

+ ^^(^78 + ^^79) + /^feo + ^e8l); 

DigiQ Dqiq Di^iQ 

= Y ^—(^32 - ^C33) + Y ^"^^^46 - ?447) + V "^^^^ ~ 

+ ^^^(64 - iim) + ]j^^{^68 - i^69) + Y^^f^(^82 + iC83) 

+ ^^(64 + ie8a) + Y^(^- + ^^-) + V^^^«^ + '^^^^ 

+ ^^^(60 + + ^^^(62 + ^63); 

+ /D^es + ^^(64 + ^e95) + ]l^iC96 + ^C97) + /^(68 + ^69); 
^15 = y^^(64 - i^s) + y^^(C66 - ^^67) + y^^(C82 " ^^83) 
+ ^j^^{i9A - ii95) + ^^15,15^6 + ^ ^^'^^ (60O + ^6oi) 

+ ^^^^(6o2 + ^6o3); 
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•7^10 — 

+ 
+ 

+ 
+ 

^■13 = 

+ 
+ 
+ 

^12 = 



(62 - «63) + y^^(C56 - ^^57) + y^^(C84 - ?C85) 



^5,14 



^2'^^ (600 - ^ClOl) + \l ^^2^'^ (^104 - ^605) + a/ ^14,14^8 

^^^Y^(6io + + y^^(6i2 + ^613); 

- + \/^(«74 - + - 

^^^y^(6o6 - ^607) + y^^(6io - + a/^13,136 

^^y^(6i4 + «6i5); 

^8,12.^ X , D1Q12. Di2,n, , 

^(^76 - 1^77) + Y ^—(^88 - ^^89) + Y (6l4 " ^Clis) 

+ ^^12,1260 + Y^ ^ 2'^^ (C1I6 + «Cll7); 

•^11 = ^^1^(66 - ^67) + - iUl) + y^^^(^78 - ^^79) 



Ad 

+ 



16,11 -t N I /-^lO'll/'t -t N I /-^14,ll/(^ V 

^(^90 - + Y (?108 - K109) + Y (?112 - K113) 

I 

w y^(6i6 - ^617) + a/-Dii,ii6i; 
^^(^80 - iCsi) + /^(68 - i69); 

y^(e92 - ^63) + ^^(602 - ^^103). (B.89) 



In numerical simulations, we have a factor ]^tAz Langevin noises and 
^^^^ for correction terms. 
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APPENDIX C 



MULTIMODE DESCRIPTION OF 



CORRELATED TWO-PHOTON STATE 



In this Appendix, we introduce a general model for quantum detection efficiency 
for multimode analysis in various quantum communication scheme. Based on this 
detection model with the spectral description of correlated two-photon state, we derive 
the effective density matrix conditioning on the detection events of entanglement 
swapping, polarization maximally entangled (PME) state projection, and quantum 
teleportation. 



To account for quantum efficiency of detector and the affect of its own spectrum 
filtering, we introduce an extra beam splitter (B.S.) with a transmissivity 77(0;, wq) 
[117] before the detection event. 77 models the quantum efficiency of the detectors 
in the microscopic level (response at frequency cuq) and the macroscopic level (time- 
integrated detection) . One example of conditioning on the single chck of the detector, 
the output density operator becomes 



C. 1 Quantum Efficiency of Detector 




1 




(C.2) 



(C.3) 
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Figure C.l: Model of quantum efficiency of detector. 



where Tr^gj is the trace over the reflected modes mg, and the flat spectrum projection 
operator IIi (only photon number is projected and no frequency resolution) is con- 
sidered in the measurement process [76]. In Figure C.l, ra\ is the incoming photon 
operator before the detection, m\ is the reflected mode, and ni\ is now the detec- 
tion mode with a modelling of spectral quantum efficiency and an effective quantum 
efficiency is defined as 



C.2 Multimode Description of Entanglement Swap- 
ping 

From Eq. (5.4), we use single mode $((X') for Raman photon and a multimode de- 
scription f{us,0Ji) for cascade photons and rewrite the effective state. Note that a 
symmetric setup is considered so the mode description is the same for both sides A 
and B in the scheme of entanglement swapping. 




(C.4) 
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|*)e// = m(l-r?2)x 



VV2{1 - 772) J f{uJs,uJi)dw,dwiXdl'^{uj,)dl'^{uJi) J ^uj')dw'dl'^{uj')S],\0)+ 
VvMl - - V2) j ^{uj)du:dl^\u:)S\ j /K,a;0at'^(a;;)aP(a;:)cia;;cia;:|O). 

(C.5) 

With the B.S., we have aj'^ = = = = 

where is the creation operator for idler photon and d\. is for Raman photon. The 
input density operator is pj„ = |^)e//(^'| and conditioning on the pair of single 
click (^«| 2) '^1,2)) we are able to generate maximally entangled singlet or triplet state 
\^)dlcz — ^ |Q)A,g- Without loss of generality, we consider a triplet state along 
with a pair of clicks (mj,ni) and use the model of quantum efficiency in Eq. (C.l) 
with tracing over the detection modes {rh\, 77,4). Note that rh\ = y/1 — r/rfig + y^7Ti4 
and n\ = \J\ — rjfil + ■sff]n'/^ as we model the quantum efficiency in the previous 
Section. 



Pout — / d'^o'^''^m4,M{'^^m3,n3[UBS^BsPinUBS^Bs]^^ (C-6) 
J —00 

M4,4 = (/L - |0)m4(0|) ® 10)^2(01 (/I4 - |0)n4(0|) (8) |0)„2(0| (C.7) 

where the unitary B.S. operator is denoted by both sides (A and B) and NRPD 
projection operators arc used [99]. These operators project the state with single 
click of the detected mode without resolving the number of photons. / is identity 
operator. The un-normalized output density operator after tracing out these modes 
becomes 
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Pout — ^ X 

P(a;,)^i|0)(0|5Aaf K)) + J J f{cus,cui)dws^*{uji)veff{c0i)d^ 



+ 



at'^(a;,)^l|0)(0|^saf K)) | + p^^, (C.8) 

where r]eff{u!) is introduced after integration of ouq, and we denote it as an effective 
quantum efficiency for idler field cui or Raman photon at frequency cu (wavelength 
780 nm for D2 line of Rb atom), p'^^^ includes the terms that won't survive after 
the interference of telecom photons in the middle B.S. (conditioning on a single click 
of detector). They involve operators hke a^aP |0) {0\df Sb, aPa^ |0) {0\SaSb and 

The normalization factor is derived by tracing over the atomic degree of freedom. 



4 

?7i?;2(l -r;i)(l -r/2) 
2 



do^Juir]eff{uJi)\f{uJs,<^i)\'^ J duj'Jujlr]eff{uj'i)\f{u'^,uj'i}\'^+ 
dusduir]eff{oJi)\f{us,oJi)\^ j \(^\^{u)r]eff{oj)dw+ 



j \^\uj)r],ff{uj)du j \^\uj')r],ff{J)dJ (C.9) 

which will be put back when we calculate the heralding and success probabilities. 

Next we interfere telecom photons with B.S. that aj'^ = ^s'^ ~ 

and again a quantum efficiency 77(0;, cuq) for telecom photon is introduced. Use 
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4 = ~ + v^^4 trace over the reflected mode £3 conditioning on the chck 
of C4 from NRPD. The effective density matrix becomes 



42) 

Pout 



/oo 
duoTr,^ { Tr,3 sP-^b? ] ^4 } 
-00 

/oo 

da;opS(^o), (CIO) 
-00 

&M^Tr,4P^{uJo)} (C.ll) 

M4 = (4 - |0)c4(0|) ® |0)e2(0|, (C.12) 



y da;,(l-7^K))/(s,i)r(s,0 1 cia;;/(s',OA^K)clK)|0)(0|x 

J du;'J{s\t')V^)dl{u;'M{0\ J rfc^r^K") vW)/*(^'", ^')+ 
J duj'^(l-rj(u'^))f(s',i')ns',i') J dujJ(s,i)y/^)cl(ujM{0\x 
J dc^;'c4K)^^K)r(«'^0 + / dw',{l-rj{u;',))f{s\t')r{s\t)x 

J dujj{s,t)^/^)ciiusm{o\ J du':c,{uj:')^/^)r{s"',t')+ 

J du',^M^)f{s',i') j da;,v/7^/(s,04(a;.)4(a;D|0)(0|x 
j du;':y^)r{s",i) I da;rv/^r(.''',i')c4(a;n£4(u;r)}+ 

mO^l^hlO^^S^ J du.VeffiuJi) J f{s,i)dUsj ns',z)dw',X 

j dujmuj)Wff{oo)^/^)c\{ujs) {sU0){0\Sb + S\\0){Q\Sa) x 
Ci{uj'^)^/r]{u'^) J J f{s,i)dus^*{uJi)rieffiuj,)duJiX 
J J /*(s^^Orfc^;$(a;Or?e//(a;Orfu;^^^K)ci(a;,)x 

[sI\0){0\Sa + Si\0){0\SB)cM)VvR)] (C.13) 
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where a brief notation for spectrum f{s, i) = f{ujs,uji) and quantum efficiency r){uj) = 
r){u!, ujq). This quantum efficiency refers to the telecom photon. We proceed to trace 
over the detected modes and the density matrix can be simphfied by interchange of 
variables in integration. 



^(2). ^ viC^-mf /" J J / / ^ / 

J dujs{l -r){uJs,uJo))fiuJs,uJi)f*{uJs,uj'i} J duj'J{uj'^,uj'i)f*{uj'^,uJi)r]{uj'^,uJo)+ 



dwsil - r]{uJs,uJo))\f{uJs,uJi)\ j (Lj'^\f{uj'^,uj'i)\ r]{uj'^,ujQ)+ 



riiri2{l - rii){l - TI2) 



8 



I j duJiTJeffiuJi) J r]{uJs,UJo)\f{uJs,UJi)\'^duJsX 

J dc^|$Hr77e//M(^;|o)(o|^5 + ^i|o)(o|^^)+ 

J J l]{uJs,UJo)f{uJs,<^i)duJs^*i<^i)Veffi'^i)duJi J f*{uJs,'^i)^i'^'i)Veffi'^'i)(^i'X 

[sU0){0\Sa + S\\0){0\Sb)^ (C.14) 

where the trace over two photon states requires the commutation relation of photon 
operators. 



Tr[mi(a;,)7nlK)|0)(0|m4K)m4K')] 

= (CIS) 

The above is the general formulation for the un-normalized density matrix condi- 
tioning on three clicks of NRPD's. We've included spectral quantum efficiency of the 
detector either for near-infrared {fjeff) or telecom wavelength {rjt = J^^r){u!,(x!o)du!o) 
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To proceed, we assume a flat and finite spectrum response {rjeffluj) — rjeff, Vti^) — 
r)t) with the range Uq G [fi — A, + A] centered at fl (near-infrared or telecom) 
and CO & [cuq — S,uJo + S]. The widths 2 A and 2S are large enough compared to 
our source bandwidth so these detection events do not give us any information of 
spectrum for our source. A perfect efficiency also means no photon loss during 
detection. Note that the integral involves multiplication of two telecom photon 
efficiency J^^r){uj,uJo)r){uj' ,uJo)dujQ — r]f{uj) that is valid if the source bandwidth is 
smaller than detector's. 

After the integration of ujq, we have 



j djjj'J{J,M)n^>i) + {^-r]t)r)t j dws\f{uj,,uj,)\' j dw',\f{uj',,uj'^\'+ 



2 



+ 



dujsf{ujs,<^i)f {uJs,<^i)>\0}{0\ + VtVeffX 

11 dwi J \f{us,uji)\^dws j da;|$(a;)|2(4|0)(0|>SB + >5l|0)(0|>SA)- 
j j f{us,Ui)dws^*{ui)dwi j f*{us,u[)^{u[)dw[ 

(4|0)(0|^^ + ^i|0)(0|^B)}. (C.16) 

C.3 Density Matrix of PME Projection and Quan- 
tum Teleportation 

In Chapter 5.4, we have the normalized density operator Poui,n the DLCZ entan- 
gled state through entanglement swapping. With another pair of DLCZ entangled 

(2) C D 

state, pout,n ■> the joint density operator for these two pairs constructs the polarization 
maximally entangled state (PME) projection and is interpreted as 
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42),AB .{2),CD _ 
Pout,n ^ Pout,n — 

(^{a^|0)(0| + ^f[|0)..(0| 

+ XiSI\0){0\Sd + XiSU0){0\Sc 
+ XiSU0){0\Sa + XiS],\0){0\Sb 
+ XiSI\0){0\Sd + XiSI\0){0\Sc 
+ XiSU0){0\Sa + XiS\\0){0\Sb 



S},\0){Q\Sc + S},\0){0\Sn 

+ \0)cd{0\[sU0){0\Sb + S\\0){0\Sa 
+ ^1(^sI\0){0\Sc + sU0){0\Sd 



SI\0){0\Sb + S\\0){0\Sa 



(C.17) 



which is used to calculate the success probability after post measurement [a click from 
each side, the side of (A or C) and (B or D)]. a = r)r{2 — rj)^l + J2j = 4, and 

Vr — Vi/V2, V — Vt, Xj is Schmidt number that is used to decompose the two-photon 
source from the cascade transition. 

In DLCZ protocol, quantum teleportation uses the similar setup in PME projec- 
tion and combines with the desired teleported state, |$) = {doS\^+ diS\j\0) , which is 
represented by two other atomic ensembles Ii and l2- The requirement of normaliza- 
tion of the state is + Mi P = 1, and the density operator of quantum teleportation 
is pQT = |$)($| ® ® pS,n^- Conditioning on clicks of Di^ and Di^, the 

effective density matrix for quantum teleportation is (using S\_^ = (Dj^ + Da)/\/2, 

= {D12 + Dc) I \/2 for the effect of beam splitter) 
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PQT,eff - 



dod^ 



(Dl^mOlDjJ + ^(Dl\0){0\D,,) + -\^{Dl^\0){0\Dj,) 



+ 



{a + by 



a'\0){0\ + 



ah 



>AB\ 



/D\\0)(0\Di, . D\ . D.\ 

^-JiLJlUl + sU0){0\Sn + Xi-^\0){0\Sn + XiSUo){0\-^) 

+ |0)cd(0| (^L|0)(0|^b + + Ai4|0)(0|^ + X,^\0){0\S, 



+ 



+ -S;|0)(0|^,, + X,-^\0){0\Sd + X^S],\0){0\^y 



. D\ |0)(0|L>/, 



+ Ai4|0)(0|^ + Ai^|0)(0|5b 



(C.18) 



which is used to calculate the success probability for teleported state. 
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APPENDIX D 



HAMILTONIAN AND EQUATION OF MOTION 
FOR FREQUENCY CONVERSION IN A 
DIAMOND TYPE ATOMIC ENSEMBLE 

In this appendix, we derive the Hamiltonian and the Maxwell-Bloch equation for 
frequency conversion in ladder- type transition. The steady state solutions for atoms 
are solved, and the solution to the field equations are discussed in Chapter 6. Similar 
to the derivation in Appendix B where the cascade emissions are investigated, the 
conversion scheme here also involves four- wave mixing with two classical driving lasers 
and two quantum fields, signal and idler. The driving lasers are applied in a way 
that signal or idler is converted only when an idler or signal is put into interaction 
with the atoms (see Figure 6.1). We will use the same quantization procedure for 
electromagnetic fields as discussed in Appendix B.2.1. 

D.l Hamiltonian and Maxwell-Bloch Equation 

To derive the coupled Maxwell-Bloch equations it is convenient to employ a quan- 
tized description of the electromagnetic field [29] and use Heisenberg-Langevin equa- 
tion methods, and then invoke a standard semiclassical factorization assumption. The 
propagation length L is discretized into 2M-I-1 elements. The positive frequency com- 
ponent of the electric field operator is given by -^+(2;) = Y^n=-M ^J~^^e^'''^'~^''"^^Cn 
where [c„, cj^,] = Snn', kn = , ujs,n = ujs + knC , n = -M, M and oug = kgC is 
the central frequency. Define the local boson operators di = ^2m+i ^n=-M CnC**^"^' 
where [di, aj,] = Sw. Similar definitions hold for the signal, s, and idler field, i, which 
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carry an additional index in the following. 

The Hamiltonian for the interacting system, Hj = —d ■ E, depicted in Figure 6.1 
is given by, (we ignore the interactions responsible for atomic spontaneous emission 
for the moment) 



H^Ho + Hi, (D.l) 

where 



3 M M 

i=l l=-M l=-M 1,1' 

M 

+ f}Wi^ aliUi^i + h^ujii>alidi^i>, (D.2) 

l=-M 1,1' 

and 



M 
l=-M 

+ gsV2M + ia'.^^ds,ie-'''''' + ft V2M + lagdi^ie'''^'^ + /i.e.} (D.3) 
where = ' '^m'" = ' ^^e Rabi frequencies ^a,{b){t) = 



fa,{b){t)dw,{2^)^{ka,(p-))/{2h) is half the standard definition, and fa,{b) is a slowly vary- 
ing temporal profile without spatial dependence (ensemble scale much less than pulse 
length). The dipole matrix element dmn = {m\d\n), coupling strength ft (j) = 
d2i,(^Q)^{ks,{i))/K £{k) = ^/h^fleoV, and Zp = 2hhl' P ^ -M,...,M. The ma- 
trix uJw = X]^-M ^"^**^"^^'~^''V(2-^ + 1) accounts for field propagation by couphng 
the local mode operators. 

The dynamical equations including dissipation due to spontaneous emission may 
be treated by standard Langevin-Heisenberg equation methods [30] , and we define 
as the natural transition rate from \ j) — )■ Since we are interested in a semiclassical 
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description, we replace the field operators by c-numbers in the Langevin equations, 
and drop the zero-mean Langevin noise sources. All atomic spin operators are also 
replaced by their expectation values. Finally, in the co-moving frame coordinates z 
and T — t — z/c the atomic equations are 



— aoi = (?Ai - ^)c^oi + i^ai^m - ^ii) + ig*sao2E~ - igia\^Ef , 

^CTo2 = (i^2 - y )5-o2 - i^i2^a + Ws^oiEt + iP*ao3Qb - iP*gi^32El, 
d 

—ail = -Toic^ii + 7i2C^22 + i^aO-Qi - i^l^oi - igsO-l^E^ + ig*ai2E~ , 
d 

^C^22 = -725-22 + igs^l2Es - ig*s^l2E~ iObO-^a - «^6C^32, 

d 

^CTgg = -7o3Cr33 + 732Cr22 " ^^^bCr^2 + ^^65-32 + igia^^E^ - ig*ao3E^ , 

^cri3 = {iAui - iAi - ^^°^ )ai3 - i^lao3 - iPgsol^E^ + iP^lai2 
+ igi^liEt, 

— ao3 = {i^i^i - ^)c^03 - i^a^i3 + iP^l^02 + igii^oo - cr33)E^, 

= {-i^b - ^°^^^^ )o-32 - iP*g*s^i3E~ + ini{a22 - ^33) + iP*gial2Et 

(D.4) 

where 72 = 712+732, P = e*^'^^"*'^'^*, the four-wave mixing mismatch wavevector AA; = 
ka — kg+kh — ki, the frequency mismatch Acu — cua+ous—oub—uji — Ai — Ab+AcUs — AuJi, 
and various detunings are defined as Aoui — uui — 003, AoUg — uig — ^^12, ^1 — ^^a ~ ^1, 
A2 = coa + ^j-'s — <-^2 = Ai + Acus , Ab = ojb — u)23- The slow-varying atomic operators 
are defined 
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1 1 1 

~ _ / -ikaZi+iuJat ~ — -I iksZi+iuJst ~ — I -ikaZi+iksZi+iUst+il^at 

~ — ^ I -iuiat+iuiit+ikaZi-ikiZi ~ — ^ I -ikiZi+iu/it ~t _ ^ l\ -iujit -ikbZi 

"^13 — C'la^ ) '-'03 — c'oac , — ^ (J^2^ e , 

_ 1 « _ 1 i _ 1 / / N 

(722 = ;^t^22> ^^33 = ^^^^33' ^11 = (D-5) 



where A^;,(2M + 1) = A^. 
The field equations are 

IS-'^^- (D.e, 

|.r^"^.o, (B.T) 

where the field operators are defined as 

E-(z,t) = V2M + lalie-'''^\ E+(z,t) = V2M + la^^^e^'^'*. (D.8) 

Langevin noises are not concerned here for we are interested in the normally- 
ordered quantity, frequency conversion efficiency, of input field and additional quan- 
tum noise corrections vanish as the |2) — )■ |3) transition driven by pump laser b has 
vanishing populations and atomic coherence. For energy and momentum conserva- 
tion {P — 1), and in the weak field limit, we solve atomic operators in steady state 
after linearizing with respect to the probe fields 



Toiaoi = il^a(l - 2aii - 022 - crsa), 

^3*2^^32 = ^^(5'22 - 5-33), 

Tis^is = -iQ*CTo3 - igs^UE'^ + i^lai2 + igial^E'l, 

7'l2Cri2 = -i^*aCFQ2 + Ws{^ll " ^22)Et + i^lZ^b-, 

7o35-o3 = -i^aCTiz + 1^02% + igi{&m - &iz)E'l (D.9) 
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where Toi ^^-iA,, T*^ = i^^H^,^ T02 ^^-1^2, ^13 = + iA,-iAuji, 

T12 = 2oi±22 jjjg = T|3 note that a"o2, 5"i3, 5"i2, 5"o3 are expressed in 

first order of fields and ctqi, a"32 in zeroth order. For population operators, we solve 
them in the zeroth order of fields and the nonzero steady states of population and 
coherence operator are (s denotes steady state solution) 



Substitute the above back into Eq. (D.9) and solve for ai2 and a"o3. The parametric 
coupling equations for the signal and idler fields become 



dz 
where 



^ E+^KiE+ + aiEt (D.ll) 



Fll,sU03 + -7^ V ) 7?, U03 + 



(D.12) 

= =— [(7oo,5(— ^ + — ^) + — ; (^03 H )J, (D.13) 

CJ-y J. 02 J- 13 -1-13 -'02 

fii = —^[<^iiA^ ^^f^) + ^ (-^12 + Tf; )J, (D.14) 

Ciy I02 -4 13 J- 02 -4 13 

«i = F00,s(Jl2 + -Tf; \- -7^) Tf. 1^12 H Tf. )l 

CJJ lQ2 J- 13 J- 13 J- 02 

(D.15) 

U = I12J-03 + -l-12[^ H — ) + -l03[^ H — ) H =^7=; . 

-'13 -^02 -4 02 -4 13 -t02-tl3 

(D.16) 
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Figure D.l: Self-coupling coefficient aj. A dimensionless quantity aiL is plotted 
with real (solid blue) and imaginery (dashed red) parts as a dependence of idler 
detuning Ao;, showing a normal dispersion inside the EIT window. 

The absorption coefficient for idler field is the real part of (— ««) and phase velocity 
Vp = oj/k{uj) = c/n{uj) that k{u) = n{u)u/c. The wavevector is related to coefficient 
tti that k{u) = Im(aj) + u/c so n{u) = 1+ lm{ai) / {u / c) . The group velocity 
is Vg = du/dk{u) = c/{n + udn/du) where n ^ 1, and it is this steep slope of 
refractive index that makes a large group delay [dn/du > inside EIT window). 
As an example in Figure D.l, we demonstrate the real and imaginary parts of self- 
coupling coefficient with the optical depth (opd) paL = 150 (see Chapter 6 for 
more details on other parameters). The dispersion curve (Im(Q;jL)) inside the left 
parametric coupling window bounded by two absorption peaks (Re(ajL)) shows a 
normal dispersion indicating a group delay at the center of the window (see Figure 
6.3 for complete parametric coupling windows). Note that we plot out unitless a^L 
where L is in the order of millimeter for regular cold atomic ensemble, and see Sec. 
II and III for detail discussion of various coupling coefficients in Eq. (D.ll) and 



162 



efficiency dependence on optical deptfi. 
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